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I.  INTRODUCTION 


Recently,  The  Technical  Cooperation  Program  (TTCP),  WTP-4  (Propulsion  Technology), 
Key  Technical  Area  10  (KTA-10)  adopted  a  lumped-parameter  interior  ballistics  model.  The 
model  is  well  documented  in  the  report  of  the  activities  of  the  KTA.1  IBRGA  was  written  in 
order  to  test  that  model,  to  compare  its  predictions  with  those  of  other  established  models,  and 
to  provide  a  relatively  simple  interior  ballistics  model  for  other  uses  by  the  community.  IBRGA 
is  a  FORTRAN  encoding  of  the  TTCP  model.  It  has  been  compiled  and  run  on  computers  as 
small  as  an  IBM  PC  and  as  large  as  a  CRAY-2.  This  report  documents  IBRGA  quite 
thoroughly,  including  comparisons  with  IBHVG2,  a  complete  FORTRAN  listing  of  the  IBRGA 
code,  and  sample  inputs  and  outputs.  A  machine-readable  copy  of  the  code,  complete  with 
sample  inputs  and  outputs,  may  be  obtained  from  the  authors  of  this  report. 

n.  DISCUSSION  OF  IBRGA 


IBRGA  performs  all  its  calculations  with  variables  expressed  in  metric  units,  specifically, 
meters,  kilograms,  and  seconds.  The  input  data  units  were  made  to  conform  to  those  of 
TDNOVA2  for  ihe  sake  of  uniformity.  Some  of  the  required  input  data  are  therefore  in 
centimeters,  grams,  and  milliseconds;  IBRGA  converts  to  the  basic  metric  units  before  the 
interior  ballistic  simulation  begins. 


There  are  two  forms  of  IBRGA  available,  the  difference  between  them  being  the  method 
used  to  calculate  the  mass  fraction  of  propellant  burned.  IBRGAS  is  a  direc)  encoding  of  the 

TTCP  model;  it  is  included  as  Appendix  A.  IBRGAS  uses  the  formula  Z,-  = 
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where  Z,  is  the  mass  fraction  burned,  v0  is  the  original  volume  of  propellant,  s(t)  is  the 

instantaneous  surface  area  and  is  the  linear  burning  rate.  IBRGAC  is  the  other  ve.sion  of 

at 

the  code;  it  uses  the  relation  Z,  =  — ,  where  Z,  is  the  mass  fraction  burned,  v(  is  the 

*0 

instantaneous  volume  of  propellant  burned,  which  is  calculated  from  the  regression  of  the 
propellant  surface  and  the  geometry  of  the  propellant  grains,  and  v0  is  the  original  volume  of 
propellant.  This  alternative  solution  technique  was  expected  to  yield  the  same  answers  but  to  be 
numerically  quicker.  The  two  versions  of  IBRGA  are  compared  in  Table  1.  Calculated  peak 
breech  pressures,  in  MPa,  are  shown  versus  the  chosen  time  step  for  the  calculations.  These 
comparisons  were  done  with  the  Cray  X-MP/48  computer,  both  for  the  long  word  length  and 
the  calculational  speed.  The  run  times  are  for  the  Cray  X-MP/48.  One  can  see  that  the  codes 
converge  to  the  same  value,  but  that  convergence  is  about  a  factor  of  100  quicker  with  IBRGAC 
than  with  IBRGAS. 

TABLE  1 


COMPARISON  OF  IBRGAC  AND  IBRGAS 


Calculated  Peak  Breech  Pressures  in  MPa 


Time  Step  (s) 

IBRGAC 

IBRGAS 

Approx.  Run  Time  (s) 

.05 

517.451 

508.526 

0.10 

.005 

517.574 

516,642 

0.44 

,0005 

517.582 

517.489 

3.90 

.00005 

517.583 

517.574 

37.40 

.000005 

517.583 

517.582 

373.00 
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Based  on  these  findings,  all  subsequent  calculations  with  IBGRA  in  this  report  used 
IBRGAC  with  a  time  step  of  0.005  seconds. 

The  basic  physical  assumptions  involved  in  IBRGAC  are  the  same  as  those  for  the  lumped- 
parameter  interior  ballistics  code  called  IBHVG23.  Both  codes  can  make  use  of  the  Lagrange 
gradient  equation  in  addition  to  other  gradients.  To  test  the  accuracy  of  IBRGAC,  comparison 
computer  runs  were  made  with  IBHVG2.  Both  codes  used  essentially  the  same  nominal  120-mm 
data  base  and  employed  the  Lagrange  gradient  equation.  Results  of  these  calculations  are 
summarized  in  Table  2.  The  percentage  difference  between  the  results  of  the  two  codes  is 
shown  to  the  right.  With  differences  between  variables  compared  of  the  order  of  0.01%,  the 
agreement  between  the  two  codes  is  quite  satisfactory. 

TABLE  2 


COMPARISON  OF  DATA  BETWEEN  IBRGAC  AND  IBHVG2 


Trajectory  variables 

IBRGAC  IBHVG2  %  difference 

Maximum  breech  pressure(MPa) 

517.574 

517.588 

0.003% 

Maximum  mean  pressure(MPa) 

460.047 

460.060 

0.003% 

Maximum  base  pressure(MPa) 

344.994 

345.003 

0.003% 

Muzzle  velocity  (m/s) 

1576.09 

1575.95 

0.009% 

Energy  variables  (J) 

IBRGAC 

IBHVG2 

%  difference 

Total  chemical  energy 

48,384,300. 

48,373,692. 

0.022% 

Energy  remaining  in  the  gas 

30,765,100. 

30,766,890 

0.006% 

Projectile  kinetic  energy 

12,173,200. 

12,164,775. 

0.069% 

Energy  lost  to  gas  and  propellant  motion 

4,059,700. 

4,056,875. 

0.070% 

Heat  energy  lost  to  convection 

1,386,220. 

1,385,152. 

0.077% 

In  addition  to  the  Lagrange  pressure  gradient  equation,  IBRGAC  permits  the  use  of  an 
alternative  gradient  equation  which  takes  into  account  the  "chambrage"  of  the  gun.  This  other 
equation  is  known  as  the  chambrage  gradient  equation  4,5  since  it  takes  into  account  the  effects 
of  a  chambrage  chamber. 

Straight  and  chambrage  chambers  are  defined  in  the  following  manner.  A  straight  or  bore- 
diameter  chamber  has  the  same  radial  distance  as  the  gun  bore.  A  chambrage  chamber  has  a 
radial  distance  which  is  larger  than  the  bore  radial  distance.  This  larger  diameter  is  tapered  or 
"necked  down",  at  the  forward  end  of  the  chamber,  to  the  same  radial  dimensions  of  the  bore 
(see  Figure  1). 

We  tested  IBRGAC  with  Lagrange  and  chambrage  gradient  equations  to  be  sure  that,  for 
data  bases  that  specified  a  bore-diameter  chamber,  a  calculation  with  the  chambrage  gradient 
gives  essentially  the  same  results  as  the  Lagrange  gradient.  Table  3  indicates  that  the 
chambrage  gradient  reduces  to  the  Lagrange  gradient  for  a  bore-diameter  chamber 
configuration. 


FIGURE  1 


A  is  the  length  of  the  tapered  section  of  the  chambrage  chamber 
B  is  the  overall  chamber  length 
C  is  the  bore  diameter 
D  is  the  chamber  diameter 
D/C  is  defined  as  the  "chambrage" 


TABLE  3 

COMPARISON  OF  CALCULATED  RESULTS  FOR  THE  LAGRANGE  AND  CHAMBRAGE 
GRADIENTS  IN  A  BORE-DIAMETER  CHAMBER 

Trajectory  IBRGAC  IBRGAC  %  difference 

(Lagrange)  (chambrage) 

Maximum  breech  pressure  (MPa)  517.58  517.62  0.01% 

Maximum  mean  pressure  (MPa)  460.05  460.08  0.01% 

Maximum  base  pressure  (MPa)  344.99  345.02  0.01% 

Muzzle  velocity  (m/s)  1576.09  1576.11  0.0012% 

III.  THE  EFFECTS  OF  CHAMBRAGE 


Table  4  demonstrates  the  effects  of  chambrage.  In  Table  4,  calculations  with  IBRGAC  are 
compared  with  those  of  XKTC  6  ,  which  is  a  phenomenologically  complete  two-phase  flow 
interior  ballistic  computer  code  which  has  a  history  of  good  agreement  with  gun  firings. 
IBRGAC  Ls  used  with  the  chambrage  gradient  equation  for  all  of  these  calculations,  since  we 
have  seen  that  it  correctly  handles  both  chambrage  and  bore-diameter  chambers.  All 
calculations  are  for  a  propellant-mass  to  projectile-mass  ratio  of  one.  The  first  line  of  the  table 
compares  calculations  for  a  nominal  120-mm  weapon  with  a  bore-diameter  chamber.  We  see 
that  the  agreement  between  IBRGAC  and  XKTC  is  quite  good.  For  the  second  line  of  the 
table,  chambrage  has  been  introduced.  The  total  chamber  volume  was  maintained,  but  7.62  cm 
of  chambrage  was  added,  so  that  the  overall  chamber  length  was  shortened  from  77.6  cm  to  54. 1 
cm.  Note  the  agreement  between  the  codes,  the  drop  in  the  maximum  breech  pressure,  and  the 
(expected)  concomitant  drop  in  the  muzzle  velocity.  The  third  line  of  the  table  results  from 
changing  the  propellant  web  so  that  XKTC  predicted  the  same  maximum  breech  pressure  as  had 
been  obtained  for  the  bore-diameter  case  in  the  first  line,  345  Mpa.  That  same  propellant  web 
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third  line  with  the  first  line,  we  see  that  a  gun  with  a  chambrage  chamber  can  achieve  a  higher 


muzzle  velocity  without  increasing  its  maximum  chamber  pressure. 
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TABLE  4 


EFFECT  OF  CHAMBRAGE  BETWEEN  XKTC  AND  IBRGAC 


XKTC 

Chamber 

chambrage 

max 

Volume 

breech 

(cm.) 

pressure 

(MPa) 

9832 

none 

345 

9832 

1.21 

310 

9832 

1.21 

345 

XKTC 

IBRGAC 

IBRGAC 

muzzle 

max 

muzzle 

velocity 

breech 

velocity 

(m/s) 

pressure 

(MPa) 

(m/s)’ 

1352 

346 

1363 

1303 

305 

1306 

1376 

347 

1390 

This  calculated  increase  in  velocity  due  to  chambrage  also  occurs  for  optimized  gun  systems, 
in  which  the  grain  dimensions  and  propellant  mass  are  varied  to  maximize  muzzle  velocity  while 
maintaining  a  given  maximum  breech  pressure.  To  illustrate,  this  point,  Appendix  D  documents 
a  calculation  in  which  a  bore-diameter  gun  was  optimized  for  velocity  for  a  maximum  breech 
pressure  of  346  MPa.  For  this  case,  a  muzzle  velocity  of  1398  m/s  was  achieved  using  8.7  kg  of 
propellant  with  a  web  of  0.2299  cm.  Appendix  E  documents  a  companion  optimization 
calculation  for  a  chambrage  chamber,  again  for  a  maximum  breech  pressure  of  346  MPa.  A 
muzzle  velocity  of  1408  m/s  was  achieved  using  8.85  kg  of  propellant  with  a  web  of  0.2276  cm. 

Many  propellant  charges  have  been  designed  in  the  past  by  using  codes  with  Lagrange 
gradient  equations  and  modeling  bore-diameter  chambers.  One  might  ask  whether  the 
propellant  dimensions  that  were  determined  from  these  calculations  were  in  fact  the  optimum 
dimensions  for  the  real,  chambrage  gun.  IBRGA  makes  this  determination  easy.  If  one  takes 
the  propellant  web  of  0.2299  cm  from  the  bore-diameter  optimization  above  and  uses  it  in  a 
chambrage  calculation,  and  further  increases  the  propellant  weight  to  8.91  kg  to  restore  346 
MPa,  one  achieves  a  muzzle  velocity  of  1408  m/s,  the  same  velocity  as  was  predicted  in  the 
optimized  chambrage  calculation  above.  Thus  we  see  that  past  optimizations  of  propellant  web 
have,  in  tact,  determined  nearly  optimal  propellant  web  for  guns  with  chambrage.  The  weight  of 
propellant  needed  was  not  so  w'ell  determined,  however. 

There  are  two  factors  inherent  in  a  chambrage  chamber  which  may  account  for  the  increased 
muzzle  velocity  which  is  both  observed  and  calculated.  First,  the  smaller  chamber  length  of  a 
chambrage  configuration  translates  into  a  smaller  pressure  gradient  for  a  significant  portion  of 
the  early  part  of  the  ballistic  cycle.  With  this  condition,  a  higher  base  pressure  and  an  increase 
in  projectile  velocity  will  result.  Second,  the  increased  radial  dimension  of  the  chamber  means  a 
larger  amount  of  expanding  gas  is  available  to  act  on  the  base  of  the  projectile,  for  a 
significantly  increased  duration,  which  would  also  lead  to  an  increase  in  energy  imparted  to  the 
projectile. 


IV.  CONCLUSIONS 

1.  IBRGA  is  a  simple  FORTRAN  computer  code  which  incorporates  the  TTCP 
lumped-parameter  interior  ballistic  model. 

2.  Calculations  of  mass  fraction  burned  which  use  the  propellant  form  function  directly 
are  more  computationally  efficient  than  those  proposed  in  the  TTCP  model. 

3.  Calculations  with  IBRGA  compare  favorably  with’ those  from  IBHVG2. 

4.  IBRGA  permits  the  use  of  a  pressure  gradient  equation  which  takes  chambrage  into 
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account;  calculations  with  IBRGA  on  the  effects  of  chambrage  compare  favorably  with 
equivalent  XKTC  calculations. 

5.  IBRGA  can  be  used  to  show  that  a  gun  chamber  with  chambrage  will,  for  the  same 
conditions,  give  a  lower  maximum  breech  pressure,  or  if  the  propellant  grain  dimensions  are 
changed  to  reproduce  the  original  maximum  breech  pressure,  will  give  a  higher  projectile 
velocity  than  a  bore-diameter  chamber. 

6.  Optimized  gun  calculations  done  with  either  the  Lagrange  gradient  or  the  chambrage 
gradient  in  a  bore-diameter  chamber,  will  give  nearly-optimal  propellant  grain  dimensions. 
Thus,  previous  optimizations  with  older  codes  led  to  proper  propellant  grain  manufacture. 
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APPENDIX  A 


Interior  Ballistic  Equations 

IBRGA  (  and  also  IBHVG2  )  relies  on  mathematical  models  in  which  only  the  most 
essential  forces  are  included.  Following  is  a  compilation  of  the  relations  which  make  up  the 
interior  ballistics  governing  equations. 

1.  The  equation  of  motion  of  the  center  of  mass  of  the  projectile  is 

,  A(Pb-br-Pg ) 


where  the  area  of  the  base  of  the  projectile  including  the  appropriate  portion  of  the  rotating 
band  is 


*  =  7d*2  * 


GLR  DG2  DL2 
GLR  +  1  +  GLR  +  1 


where  Db2  — 


The  pressure  on  the  base  of  the  projectile,  Pb,  including  the  approximate  pressure  gradient  effect  is, 


for  the  Lagrange  gradient, 


p  +  £r(PLt£A\ 

f  ,  .  CT  } 


and  for  the  chambrage  gradient 


P  -  6 


where  5=  -a\J\  -  bJ2  + 


+^4  b 

V{zP) 


and  7=1-  a2J 5  + 


V(z>)  ’ 


V(zp)  is  the  volume  up  to  the  base  of  the  projectile  , 

?,id  zp  is  the  distance  from  the  breech  to  the  base  of  the  projectile. 

CtA  Avp  2  A  ( br  +  Pg) 

0l~  V2  (zp)  V(zp)  +  mp 


-CtA2 

>npv2(zp) 
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Also, 


c-i 

11 

-0 

(110) 

1  Jo  A(z)  ’ 

(Ml) 

j, 

2  A*(z)' 

(1.12) 

J>  =  g  A(z)h(z)d2, 

(1.13) 

h  =  f*p  A(z)J2(z)dz  , 

(1.14) 

where  A  (z)  is  the  area  at  distance  z  from  the  breech, 
V(z)  is  the  volume  at  distance  z  from  the  breech, 
and  J\  ,  J2  ,  Ji  ,  J 4  are  numerically  evaluated  at  zp. 


2.  The  velocity  of  the  center  of  mass  of  the  projectile  is 


3.  The  breech  pressure  Pa  is  , 
for  the  Lagrange  gradient , 

K  =  p,  +  (3.1) 

and  for  the  chambrage  gradient 


P 0  =W  -  atJy)  -  a\J\  -  bJ 2 

4.  The  travel  of  the  projectile  is 

*  =  (vpdt  +  (v'Pdl  ■ 

5.  The  mass  fraction  burning  rate  of  the  ith  propellant  is 

Si  r, 


Z,  = 


(3.2) 


(4.1) 


(5.1) 


St 


where  r,  and  S ,  are  the  instantaneous  values  of  the  burning  rate  and  surface  area,  respectively, 
and  Vgi  is  the  initial  grain  volume. 


The  linear  burning  rate  r4  is  given  by 


(5.2) 
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6.  The  fraction  of  mass  burned  of  the  ith  propellant  is 

Z,  =  j‘  Zt  dt , 

or  Z{  =  TT~  ' 

where  v,  is  the  instantaneous  volume  of  the  propellant  burned  calculated  from 
the  distance  burned  into  the  propellant  given  by  f  rt  dt . 

7.  The  space- mean  pressure  P  (  Noble- Able  Law  )  is 

T  ["  FjCjZj  F,C, 

‘  v  2j  T’  'r  f  ’ 

c  I  «  1  01  1  01 

where  the  number  of  propellants  is  n  and  the  volume  available  for  gases  is 

K  =  V0  +  Ax  -  £  %  1  -  Zi)  -  XQbtZ,  -  C,bt , 
i  Pi  i 

and  the  temperature  of  the  gases  is  given  by 


H  F'CiZi 

F>c,  _£  _ 

F  —  F  —  F>  —  F  F  F 

r  (7.--1) 

+ (-» - 1)  " 

Epr  Zbr  Zr 

f*  FiCiZi 

F,C, 

(?  (7/  -  1  )T0I 

(7/-  1  )T0, 

The  energy  loss  due  to  projectile  translation  is 

r-  mp  VP  1 

“p1  ~  2  ’ 

the  energy  loss  due  to  projectile  rotation  is 

^  *  2  mpvp2Tw2 

EPr  ~  4 

the  energy  loss  due  to  propellant  gas  and  unburned  propellant  motion  is  , 
for  the  Lagrange  gradient , 

CrvJ 


and  for  the  chambrage  gradient 


CrVpz  A2Jq 


p  2  v\Zp)  ' 

The  energy  loss  for  work  against  bore  resistance  due  to  engraving  and  friction  Ls 
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the  energy  loss  due  to  recoil  is 


Ebr  -  A  fo  brvpdt , 


P  m'p  <yrpf 
E-  -  2  • 


and  the  energy  lost  due  to  air  resistance  is 

Ed  =  A  /Q  vpPgdt . 

The  energy  lost  due  to  heat  transfer  to  the  chamber  walls  is  given  by 


with 


where 


Q=Awh(T-T,), 


and 


Aw  =  —  xDb  +  2A  +  xDb  x  , 
A 


h  =  A  Cppv  +  h0  , 

where  the  Nordheim  heat  transfer  coefficient  A  is 

A  =  [  13.2  +  4  log  10  [  1CK).  Db  ]]-2 » 
and 

_  1 

v  -  2  *p  . 


SQZ,-  +  C/ 


P  ~ 


K 


F,  l, 


^  Cp‘  =  ( 7/  -  1 ) 


"  FjjjCjZ^  FrnQ 

r  (7/  -  i)Tei  T  (7/  -  1)7; 


01 


T.CiZi  +  C, 


(7.8) 

(7.9) 

(7.10) 

(7.11) 

(7.12) 

(7.13) 

(7.14) 

(7.15) 

(7.16) 

(7.17) 

(7.18) 

(7.19) 
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(7.20) 


Thus  the  temperature  of  the  chamber  wall  is 

+  f^br 


rp  _  ,  rp 

e~  Cpt'fiW  Aw  Dw  0  ’ 


where  T0  is  the  initial  temperature  of  the  wall. 

8.  The  equation  of  motion  for  the  recoiling  parts  (See  figure  A.)  is 

A 


where 


When 


V'P  ~ 


m 


rp 


(Po-^-br), 


vrB  =  0  until  P0  > 


RPn 


vrp  —  \j  umu  i  g  "  ^  br . 


RP 

Po  =  ~~7~  »  1  —  (r0  and  ir  =  t  -  tTo  , 


and  if 


vrp  <  0 , 


vrp  =  0  . 


9.  The  velocity  of  the  recoiling  parts  in  the  ground  reference  frame  is 


(8.1) 


(8.2) 


(8.3) 


(8.4) 


(9.1) 


WEAPON  RECOIL 


TIME  OF  RECOIL 


LIST  OF  SYMBOLS 


symbol  definition 


units 


A 

Aw 

A(z) 

bi 

bi 

br 

Q 

Q 

Cp 

CDI 


Cp  H> 

Ct 

D„ 

DG 

DL 

Dw 


Ed 

Eh 

Ep 

Ep  r 

EP' 

Er 

f 

Fi 

F, 

GLR 

h 


mrp 


area  of  base  of  projectile  including  appropriate  portion 
of  rotating  band 

chamber  wall  area  plus  area  of  gun  tube  wall  exposed  to 
propellant  gases 

area  at  a  distance  z  from  the  breech 
covolume  of  ith  propellant 
covolume  of  igniter 

bore  resistance  due  to  friction  and  engraving 
initial  mass  of  ith  propellant 
initial  mass  of  igniter . 

specific  heat  at  constant  pressure  of  propellant  gas 

specific  heat  at  constant  pressure  of  ith 

propellant  (  over  temperature  range  T  to  Tol  ) 

heat  capacity  of  steel  of  chamber  wall 

total  mass  of  propellants  and  igniter 

diameter  of  bore 

diameter  of  grooves 

diameter  of  lands 

chamber  wall  thickness 

energy  lost  to  work  against  bore  resistance  due  to 
friction  and  engraving 
energy  lost  to  air  resistance 
energy  loss  from  heat  transfer  to  chamber  wall 
energy  lost  to  propellant  gas  and  unburned  propellant 
motion 

energy  loss  due  to  projectile  rotation 
energy  loss  due  to  projectile  translation 
energy  loss  due  to  recoil 

fraction  of  work  done  against  bore  friction  that  preheats 
chamber 

force  per  unit  mass  of  ith  propellant 
force  per  unit  mass  of  igniter  propellant 
groove  to  land  ratio 

heat  transfer  coefficient  of  Nordheim,  Soodak,  and 

Nordheim 

mass  of  projectile 

mass  of  recoiling  parts 


m 

m3/kg 

m3/kg 

Pa 

kg 

kg 

J/kg-K 

J/kg-K 

J/kg-K 

kg 

m 

m 

m 

m 


J 

J 

J 

J 

J 

J 

J 


none 

J/kg 

J/kg 

none 


2 

watt/m 

kg 

kg 


; 


I 


a 


a 


l 


I 


i 
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1 


P  space  mean  pressure 

Pb  pressure  on  base  of  projectile 

Pg  pressure  of  gas  or  air  ahead  of  projectile 

P0  breech  pressure 

Q  heat  flux  to  the  chamber  wall 

r,-  linear  burning  rate  of  ith  propellant 

RP  resistive  force  to  projectile  motion 

Sf  surface  area  of  partially  burned  ith  propellant  grain 

t  time 

tr  recoil  time 

T  mean  temperature  of  propellant  gases 

Tc  temperature  of  chamber  wall 

T0i  adiabatic  flame  temperature  of  ith  propellant 

TOJ  adiabatic  flame  temperature  of  igniter  propellant 

T0  initial  temperature  of  chamber  wall 

Tw  twist  of  rifling 

V(z)  volume  at  distance  z  from  the  breech 

V(zp)  volume  up  to  the  base  of  the  projectile 

v  mean  gas  velocity 

vp  velocity  of  projectile 

vp  acceleration  of  projectile 

vrp  velocity  of  recoiling  parts 

vc  volume  behind  projectile  available  for  propellant  gas 

volume  of  an  unburned  ith  propellant  grain 
va  volume  of  an  empty  cannon  chamber 

v,  instantaneous  volume  of  propellant  burned 

x  travel  of  projectile 

2/  fraction  of  mass  burned  of  the  ith  propellant 

Z,  mass  fraction  burning  rate  for  ith  propellant 

zp  distance  from  the  breech  to  the  base  of  the  projectile 

a,  burning  rate  exponent  for  the  ith  propellant 

Pi  burning  rate  coefficient  for  the  ith  propellant 

7 i  ratio  of  specific  heats  for  ith  propellant 

7/  ratio  of  specific  heats  for  igniter 

A  Nordheim  heat  transfer  coefficient 

p  mean  gas  density 

Pi  density  of  ith  propellant 

pw  density  of  chamber  wall  steel 


Pa 

Pa 

Pa 

Pa 

w 

m/s 

N 

m2 

m 

s 

s 

K 

K 

K 

K 

K 

turns/caliber 

3 

m 

3 

m 

m/s 

m/s 

m/s2 

m/s 


m 

none 

s"1 

m 

none 

m/s  -Pa*1 

none 

none 

none 

kg/rr3 

kg/m3 

kg/m3 
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APPENDIX  B 


Source  code  listing  of  the  IBRGAC  program, 
program  ibrgac 

character  bdfile*10/outfil*10 

dimension  br(10) ,trav(10)  ,rp(10) ,tr(10) , forcp(lO) ,tempp(10) ,covp( 
&10) 

dimension  chwp(10) ,rhop(10) ,gamap(10) ,nperfs(10) ,glenp(10) ,pdpi(lO 
&)  ,p  dpo(10) ,gdiap(10) ,dbpcp(10) ,alpha(10,10) ,beta(10,10) , 

&pres(10, 10) 

dimension  a(4) ,b(4) ,ak(4) ,d(20) ,y(20) ,p(20) ,z(20) ,frac(10) , surf (10 
&) ,nbr(10) ,ibo(10) 
real  lambda, jlzp, j2zp, j3zp, j4zp 
dimension  chdist(5) ,chdiam(5) , bint (4) 
call  gettim(ihr,imin,isec,ihuns) 
pi=3. 14159 
write(*, 15) 

format ( 7  input  name  of  data  file  to  be  used  as  input  ') 
read ( * , 10) bdf ile 
format (alO) 

open(unit=2, err-999, file=bdf ile, status=*'old' , iostat=ios) 
rewind  2 
write (*,25) 

format ( '  input  name  of  output  file  ' ) 
read ( * , 10) outf il 

open (unit=3 , err=998 , f ile=outf il , status='new' ) 

write (3, 16) bdf ile 

format('  THE  INPUT  FILE  IS  ',alO) 

read ( 2 , * , end~2  0 , err=3  0 ) cham , grve , aland , glr , twst , travp , igrad 
if (igrad.gt. l)go  to  51 
write (3, 55) 

format (lx, 'using  Lagrange  pressure  gradient') 
go  to  52 

define  chambrage  assumes  nchpts=number  of  points  to  define 
chamber  >  or  =  2  <  or  =  5  (?) ,chdiam(I)  defines  chamber  diameter 
at  chdist  (I)  chamber  distance,  chdiam(nchpts)  is  assumed  to  be 
the  bore  diameter  and  chdist (i)  is  assumed  to  be  0,  i.e.  at  the 
breech.  Assumes  truncated  cones, 
write (3,47, err=3Q) 

format (lx, 'Using  chambrage  pressure  gradient') 
read(2,*,end=20,eri-30)nchpts, (chdist(I) ,chdiam(I) ,I=l,nchpts) 
write(3,53,err=30) (chdist(I) ,chdiam(I) ,I=l,nchpts) 
format (///,'  chamber  distance  cm  chamber  diameter  cm' ,/ (5x,el4 
& . 6 , 5x, el4 . 6) ) 
do  54  I=l,nchpts 
chdist (I) =0.01*chdist(I) 
chdiam(I)=0.01*chdiam(I) 
calculate  chamber  integrals  and  volume 
if (nchpts.gt.5)  write(3, 44,err=30) 
format (lx, 'use  first  5  points') 
if (nchpts.gt.5) nchpts=5 
bore=chd iam ( nchpt s ) 

if (chdist (1) .ne. 0.0) write (3, 45,err=30) 
format(lx,'  #  points  ?  ') 


chdist(l)«0.0 
pi3-pi/3 . 0 
bl-0.0 
b2-0.0 

b:  -o.o 

b4*0.0 

points-25.0 

56  points-points+points 

step-chd is t ( nchpts ) /points 

2Z“0.0 

bint (1) -0.0 

bint(3)-0.0 

bint (4) =0.0 

bvol=0 . 0 

r  2 =0 . 5  *chd  iain  ( 1 ) 

k=l 

j « int ( po ints+0 . 5 ) 
do  57  1=1 ,j 
zz=zz+step 

if (k.eq.nchpts-l)go  to  46 
do  58  I l=k, nchpts- 1 

if (zz.gt.chdist(Il) .and.  zz.lt.chdist(ll+l) )go  to  59 

58  continue 
Il=nchpts-1 

59  k=Il 

46  diam=(zz-chdist(k)  )/(chdist(k+l)-chdist(k) ) 

diam=chdiam(k)+diara*  (chdiam(k+l)  -chdiam(k) ) 
rl=0.5*diam 

area=pi*  ( r  l+r2 )  *  ( r  l+r!j )  /4 . 
bvol=bvol+step*pi3  * ( rl*r l+rl*r2+r2  *r2 ) 
bint ( 1 ) =bint ( 1 ) +step*bvol/area 
bint ( 3 ) =bint ( 3 ) +step*area*bint (1) 
bint (4 ) =bint (4 ) +step*bvol*bvol/area 
57  r2=rl 

temp=abs ( 1 . 0-bl/bint ( 1 ) ) 

if (abs ( 1 . 0-b3/bint ( 3 ) ) . gt . temp) temp=ahs ( 1 . 0-b3/bint ( 3 ) ) 

if ( abs ( 1 . 0-b4/bint  ( 4 ) ) . gt . temp) temp=abs ( 1 . 0-b4/bint ( 4 ) ) 

if (temp. le. 0.001) go  to  41 

bl=bint(l) 

b3=bint(3) 

b4=bint(4) 

go  to  56 

4 1  cham=bvol * 1 .  e6 

c  write(3 ,47,err=30)bint (1) ,bint(3) ,bint(4) 

c  format (lx, 'bint  1  =  ',el4.6,'  bint  3  =  ',el4.6,'  bint  4  =  ',el4. 

C  &6) 

chml en=chd is t ( nchpts ) 

52  write (3,40, err=3  0) cham , grve , aland , gl r , twc t , travp 

40  format (lx, 'chamber  volume  cm**3' ,el4 .6,/'  groove  diam  cm',el4.6,/ 

&'  land  diam  cm',el4.6,/'  groove/land  ratio' ,el4 . 6,/'  twist  turns 
&/caliber  ',el4.6,/'  projectile  travel  cm' ,el4.6///) 
cham-'Cham*  1 .  e-6 
grve=grve*l. e-2 
aland=aland*l . e-2 
travp=travp* 1 . e-2 
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read ( 2 , * , end*20 , err»30 )  prwt , iair ,  htf r , pgas 
write ( 3 , 50 , err=3  0) prwt , iair , ht f r , pgas 
50  format ( lx, 'projectile  mass  kg' , 614.6,/'  switch  to  calculate  energ 

&y  lost  to  air  resistance  J',i2,/'  fraction  of  work  against  bore  u 
&sed  to  heat  the  tube' ,el4.6/lx, '  gas  pressure  Pa'  ,el4.6) 

read ( 2 , * , end=2  0 , err=30 ) npts  , (br(i) ,trav(i) ,i=l,npts) 
write(3, 60,err=30)npts, (br(i) , trav(i) , i=l,npts) 

60  format (lx, 'number  barrel  resistance  points' , i2,/'  bore  resistance 
&  MPa  -  travel  an'/(lx,el4.6,el4.6) ) 
write(3, 65) 
do  62  i=l,npts 
br (i)=br(i) *l.e6 
trav ( i) =trav ( i ) * 1 . a-2 
62  continue 

65  format (lx) 

read (2 , * , end=20 , err*30) rcwt , nrp, (rp(i) ,tr(i) ,i=l,nrp) 
write(3,70,err=30)rcwt,nrp, (rp(i) ,tr(i) ,i=l,nrp) 

70  format (lx, '  mass  of  recoiling  parts  kg',el4.6,/'  number  of  recoi 

&1  point  pairs', i2,/'  recoil  force  N','  recoil  time  sec'/, (lx,el4 
& . 6 , 3x,el4 . 6) ) 
write(3, 65) 

read ( 2 , * , end=2  0 , err=3  0 ) ho , tshl , cshl , twal , hi , rhocs 
write  (3,75,  err=3  0 )  ho ,  tshl ,  cs. .  ll ,  twa  1 ,  hi ,  rhocs 
75  format (lx, '  free  convective  heat  transfer  coefficient  w/cm**2  K', 

&el4.6,/'  chamber  wall  thickness  cm',el4.6,/'  heat  capacity  of  st 
&eel  of  chamber  wall  J/g  K',el4.6,/'  initial  temperature  of  chambe 
&r  wall  K',el4.6,/'  heat  loss  coefficient' ,el4. 6,/'  density  of  ch 
Samber  wall  steel  g/cm**3 ' ,el4 . 6//) 
ho=ho/l.e-4 
tshl=tshl*l.e-2 
cshl=cshl*l . e+3 
rhocs=rhocs*l . e-3/1 . e-6 

read (2 , * , end=20, err=30) forcig, covi , tempi , chwi ,gamai 
write (3,85, err=30) forcig , covi , tempi , chwi , gamai 
85  format (lx, '  impetus  of  igniter  propellant  J/g' ,el4. 6,/'  covolume 
&  of  igniter  cm**3/g' ,el4 .6,/'  adiabatic  flame  temperature  of  igni 
&ter  propellant  K',el4.6,/'  initial  mass  of  igniter  kg',el4.6,/'  r 
&atio  of  specific  heats  for  igniter' ,el4. 6//) 
forcig=forcig*l.e+3 
covi=covi*l. e-6/l.e-3 

read ( 2 , * , end=2  0 , er r=3  0 ) nprop , (forcp(i) ,tempp(i) ,covp(i) ,chwp(i) , 
&rhop(i) ,gamap(i) ,nperfs(i) ,glenp(i) ,pdpi(i) ,pdpo(i) ,gdiap(i) ,dbpcp 
&(i) , i=l, nprop) 

write(3,95,err= 30) (i,forcp(i) ,tempp(i) ,covp(i) ,chwp(i) 

&,rhop(i) ,gamap(i) ,nperfs(i) ,glenp(i) ,pdpi(i) ,pdpo(i) ,gdiap(i) ,dbpc 
&p(i) ,1=1, nprop) 

95  format(('  for  propellant  number' , i2,/'  impetus  of  propellant  J/g 
&' ,el4.6,/'  adiabatic  temperature  of  propellant  K',el4.6,/'  covol 
&ume  of  propellant  cm**3/g' ,el4.6/'  initial  mass  of  propellant  kg' 
&,el4.6/'  density  of  propellant  g/cm**3' ,el4.6/'  ratio  of  specif i 
&c  heats  for  propellant' , el4 . 6/ '  number  of  perforations  of  propell 
&ant',i2/'  length  of  propellant  grain  cm',el4.6/'  diameter  of  inn 
&er  perforation  in  propellant  grains  cm',el4.6/'  diameter  of  outer 
^perforation  of  propellant  grains  cm',el4.6/'  outside  diameter  of 
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Spropellant  grain  cm',el4.6/'  distance  between  perf  centers  cm', el 
&4.6)//) 
tmpi*»0 . 0 
do  96  i»l,nprop 
f orcp ( i ) *f orcp ( i ) * 1 . e+3 
covp ( i ) *covp ( i ) *1 . e-6/1 . e-3 
rhop ( i ) «rhop ( i ) * l . e-3/ 1 . e-6 
glenp ( i) =glenp ( i) *0 . 01 
pdpi(i)-pdpi(i)*0.01 
pdpo  (i)  =pdpo  (i)  *0 . 01 
gdiap ( i) *gdiap ( i) *0 . 01 
dbpcp ( i ) =dbpcp ( i ) *0 . 01 
tmpi=tmpi+chwp ( i) 

96  continue 
tmpi=tmpi+chwi 
do  97  j=l,nprop 

read(2/*fend=20,err-30)nbr(j) , (alpha(j,i) ,beta(j ,i) ,pres(j,i) , 
&i=l,nbr(j) ) 

write(3, 110,err>*30)nbr  (j) , (alpha (j , i) ,beta( j , i) ,pres( j , i) , 
&i=l,nbr(j) ) 

110  format (lx, ' number  of  burning  rate  points' , i2/3x, '  exponent' ,8x, ' 

&  coefficient' , lOx, '  pressure'/5x, '-' , I5x, 'ciVsec-MPa**ai' , lOx, 'MP 
&a' , / (lx, el4 .6, 5x,el4 .6, 15x,el4 .6) ) 
do  112  i*l,nbr(j) 
beta(j ,i)»beta(j ,i)*l.e-2 
pres(j , i)-pres(j , i) *l.e6 

112  continue 

97  continue 
write(3,65) 

read ( 2 , * , end=2 0 , err* 30 ) del tat , del tap , tstop 
write (3,120, err*3  0 ) del tat , del tap , tstop 

120  format (lx, 'time  increment  msec' ,el4.6, '  print  increment  msec',el4 
&. 6/lx, 'time  to  stop  calculation  msec  ',el4.6) 
write(*, 130) 
deltat=deltat*0 . 001 
deltap=deltap*0 . 001 
tstop=tstop* .001 

130  format (lx, 'the  data  has  been  read') 
if (igrad.gt. l)go  to  131 

bore= (glr*grve*grve+aland*aland) / (glr+l . ) 
bore=sqrt ( bore ) 

131  areab=pi*bore*bore/4 . 

lambda=l./( (13.2+4.*logl0(i00. *bore) )**2) 

pmaxro=0 . 0 

pmaxbr=0 . 0 

pitiaxba=0 . 0 

tpmaxm=0 . 0 

tpmaxbr=0 . 0 

tpmaxba=0 . 0 

tpmax=0 . 0 

a(l)-0.5 

a ( 2 ) =1 . -sqrt ( 2 . ) /2 . 
a(3)=l.+sqrt(2.)/2. 
a(4)=l./6. 
b ( 1) =2  . 
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b(2)=l. 
b(3)«l. 
b ( 4 ) =2 . 
ak(l)*0. 5 
ak(2)=a(2) 
ak(3)=a(3) 
ak ( 4 ) =0 . 5 
vp0=0 . 0 
tr0=0 . 0 
tcw=0.0 

do  5  i=l,nprop 
ibo(i)*=0 

5  vpO=chwp(i)/rhop(i)+vpO 
volgi=cham-vpO-chw i*covi 
pmean=  fore ig*chwi/ vo lg i 
volg=*volgi 
volgi=volgi+vp0 
wallt==twal 
ptime=0. 0 
ibrp=8 
z(3)=l. 

nde= ibrp+nprop 

wr ite ( 3 , 13 2 ) areab , pmean , vpO , volgi 

132  format (lx, 'area  bore  mA2  ',el6.6,'  pressure  from  ign  Pa',el6.6,/ 
&lx, '  volume  of  unbumt  prop  mA3  ',el6.6,'  init  cham  vol-cov  ign  m 
&A3  ' ,el6.6) 
write (3, 6) 

6  format  (lx, '  time  acc  vel  dis  mpress 

&  phase  pbrch  ' ) 

iswl=0 

19  continue 

do  11  J=l, 4 

c  FIND  BARREL  RESISTANCE 

do  201  k=2,npts 

if (y(2)+y(7) .ge.trav(k) )go  to  201 
go  to  203 

201  continue 

k=npts 

203  resp=(trav(k)-y(2) -y (7) )/(trav(k) -trav(k-l) ) 
resp=br(k) -resp*(br(k) -br(k-l) ) 

C  FIND  MASS  FRACTION  BURNING  RATE 

do  211  k=l,nprop 
if (ibo(k) .eq. l)goto211 

call  prf017(pdpo(k) ,pdpi(k) ,gdiap(k) ,dbpcp(k) ,glenp(k) ,surf (k) ,fra 
&c(k)  , y(ibrp+k)  ,nperfs(k)  ,u) 
if (surf (k) .lt.l.e-10)  ibo(k)=l 

211  continue 

k=nprop 

c  ENERGY  LOSS  TO  PROJECTILE  TRANSLATION 

elpt=prwt*y(l) *y (l)/2 . 

c  ENERGY  LOSS  DUE  TO  PROJECTILE  ROTATION  . 
elpr=pi*pi*prwt*y ( 1) *y ( 1) /4 . *twst*twst 

c  ENERGY  LOSS  DUE  TO  GAS  AND  PROPEIJANT  MOTION 

if (igrad. le. l)go  to  214 
pt=y ( 2 ) +y ( 7 ) 
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vzp=bvol+areab*pt 

j4zp=bint (4) +  ( (bvol+areab*pt)  **3-bvol**3  )/3  ./areab/areab 
elgpm“tinpi*y  (1)  *y(l)  *areab*areab*j42p/2./vzp/vzp/vzp 
go  to  216 

214  elgpm=tmpi*y(l)*y(l)/6. 
c  ENERGY  LOSS  FROM  BORE  RESISTANCE 
21t  elbr=y(4) 

z ( 4 ) =areab*resp*y ( 1 ) 

C  ENERGY  LOSS  DUE  TO  RECOIL 
elrc=rcwt*y(6) *y{6)/2. 

C  ENERGY  LOSS  DUE  TO  HEAT  LOSS 

areaw*=cham/areab*pi*bore+2 . *areatH-pi*bore* (y (2)+y (7) ) 

avden*0 . 0 

avc=*0.0 

avcp*0 . 0 

zl8=0 

zl9=0 

do  213  k=l , nprop 

zl8=forcp(k) *garoap(k) *chwp(k) *frac(k)/ (gamap(k)-l. )/tempp(k)+zl8 
zl9=chwp (k) *frac (k) +zl9 
avden=avden+chwp (k) *frac (k) 

213  continue 

avcp=(zl8+forcig*gamai*chwi/ (gamai-1. ) /tempi)/ (zl9+chwi) 

avden*=  (avden+chwi)  /  (volg+covl) 

awel*.5*y(l) 

htns=l  ambda*avcp*avden*awel+ho 
z ( 5) =areaw*htns* (tgas-wallt) *hl 
elht*y(5) 

wallt= (elht+htfr*elbr) /cshl/rhocs/areaw/tshl+twal 
c  wr  it  e  ( 3 ,  * )  lambda ,  avcp ,  avden ,  awel ,  ho ,  areaw ,  htns ,  tgas ,  wal  1 1 ,  hi ,  z  ( 5 ) 

c  & , elht 

c  ENERGY  LOSS  DUE  TO  AIR  RESISTANCE 

air=iair 

z ( 8 ) =y ( 1 ) *pgas*air 
elar=areab*y(8) 
c  RECOIL 

z (6) =0. 0 

if (pbrch.le.rp(l)/areab)go  to  221 
rfor=rp(2) 

if (y (3) -tr0.ge.tr(2) )go  to  222 
rfor=(tr(2)-(y(3)-tro) )/ (tr(2) -tr(l) ) 
rfor=rp(2)-rfor*(rp(2)~rp(l) ) 

222  z ( 6) =areab/rcwt* (pbrch-rfor/areab-resp) 
if (y( 6) .It. 0.0) y( 6) “0.0 

z  (7)  =y  (6) 
goto  223 
221  tr0=y ( 3 ) 

223  continue 

c  CALCULATE  GAS  TEMPERATURE 
eprop=0. 0 
rprop=0 . 0 
do  231  k=l, nprop 

eprop=eprop+forcp('k)  *chwp(k)  *frac  (k)/ (gamap(k) -1 . ) 
rprop=rprop+forcp(k) *chwp(k) *frac(k)/(gamap(k) -1 . ) /tempp(k) 

231  continue 
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tenergy=elpt+elpr+elgpm+elbr+elrc+elht+elar 

tgas= (eprop+f orcig*chwi/ (gamai-l. ) -olpt-elpr-elgpm-elbr-elrc-elht 
fi-elar) / ( rprop+f orcig*chwi/ (gamai-l. ) /tempi) 
c  FIND  FREE  VOLUME 

Vl=0.0 
COV 1=0.0 

do  241  k=l,nprop 

vl=chwp (k) * ( 1 . -f rac (k) ) /rhop (k) +vl 
covl=covl+chwp(k) *covp(k) *frac(k) 

241  continue 

volg=volgi+areab*  (y  (2)+y(7) )  -vl-covl 
C  CALCULATE  MEAN  PRESSURE 

r 1=0.0 

do  251  k=l,nprop 

rl=rl+forcp(k) *chwp(k) *frac(k)/terapp(k) 

251  continue 

pmean-tgas/volg* ( rl+  f or cig*chwi/ tempi ) 

259  resp=resp+pgas*air 

if (igrad.le.l)go  to  252 
if (iswl.ne.O)go  to  253 
pbase=pmean 
pbrch=pmean 

if (pbase.gt.resp+l. ) iswl=l 
go  to  257 

c  USE  CHAMBRAGE  PRESSURE  GRADIENT  EQUATION 

253  j lzp=bint (1)  + (bvol*pt+areab/2 . *pt*pt) /areab 
j  2zp= (bvol+areab*pt) **2/areab/areab 

j  3zp=bint (3 ) +areab*bint ( 1 ) *pt+bvol*pt*pt/2 . +areab*pt*pt*pt/6 . 

a2t=-tmpi*areab*areab/prvt/vzp/vzp 

alf=l.-a2t*jlzp 

alt=tmpi*areab* (areab«v (l) *y (l)/vzp+areab*resp/prwt)/vzp/vzp 
bt=-tmpi*y(l) *y (1) *areab*areab/2 ./vzp/vzp/vzp 
bata=-alt*j lzp-bt*j2zp 
gaimna=alf+a2t*  j  3zp/vzp 
delta=bata+alt*j3zp/vzp+bt*j4zp/vzp 
c  calculate  b*se  pressure 

pbase= (pmean-dalta) /gamma 
c  calculate  breech  pressure 

pbrch=alf*pbase+bata 
go  to  254 

C  USE  LAGRANGE  PRESSURE  GRADIENT  EQUATION 

252  if (iswl.ne. 0)go  to  256 

if (pmean. It . resp) resp=pmean 
c  CALCULATE  BASE  PRESSURE 

256  pbase=(pmean+tnpi*resp/3 ./prvt) / ( 1 . +tmpi/3 ./prwt) 
if (pbase.gt. resp+1. ) iswl=l 

C  CALCULATE  BREECH  PE S SURE 

pbrch=pbase+tmpi/2 . /prwt* (pbase-resp) 
c  CALCULATE  PROJECTILE  ACCELERATION 

254  z (l)=areab* (pbase-resp) /prwt 
if (z(l) .lt.0.0)go  to  257 

go  to  258 

257  if  (iswl.eq.  0)  z-(l)=0,0 

258  if  (y(l) . lt.o.o)y(l)=o.o 
z(2)=y(l) 
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C  GET  BURNING  RATE 
do  264  m=l,nprop 
z(ibrp+n)  =0.0 
if (ibo(m) .eq. 1)  goto  264 
do  262  k«l,nbr(m) 
if  (pmean.gt.pres(in,k)  )go  to  262 
go  to  263 

262  continue 
k=nbr(m) 

263  2  ( ibrptm)  -beta  (m,  k)  *  (pmean*l .  e-6)  **alpha  (m,  k) 

264  continue 

do  21  i=l,nde 
d(i)=(z(i)-b(j)*p(i))*a(j) 
y ( i) =deltat*d( i) +y ( i) 
p(i)=3.*d(i)-ak(j)*z(i)+p(i) 

21  continue 
11  continue 
t=t+deltat 

if(pmaxm.gt.pmean)goto  281 

pmaxm=pmean 

tpmaxm=y(3) 

281  if (pmaxba.gt.pbase)go  to  282 
pmaxba=pbase 
tpmaxba=y(3) 

282  if (pmaxbr.gt.pbrch)go  to  283 
pmaxbr=pbrch 
tpmaxbr=y(3) 

283  continue 

if (y(3) .lt.ptime)go  to  272 
ptime=ptime+deltap 

write (3, 7) y  (3)  ,z(l)  ,y (1)  ,y(2)  ’maean , pbase , pbrch 

7  format(lx,7ell.4) 

316  format (lx) 

272  continue 

if (t.gt.t stop) goto  200 

if (y(2) .gt.travp)go  to  200 

rravelo=y(l) 

tmvelo=y(3) 

disto=y (2) 

go  to  19 

200  write(3,311)t,y(3) 

311  format(lx,'  deltat  t',  el4.6,  '  intg  t',el4.6) 
write  (3,312)  pmaxm ,  tpinaxm 

312  format ( lx, 'FWAXMEAN  Pa  ',614.6, '  time  at  PMAXMEAN  sec  ',el4.6) 
write (3,313) pmaxba , tpmaxba 

313  format ( lx , ' FMAXBASE  Pa  ',el4.6,'  time  at  PMAXBASE  sec  ',el4.6) 
write (3,314) pmaxbr , tpmaxbr 

314  format ( lx, 'FMAXBREECH  Pa  ',el4.6,'  time  at  PMAX BREECH  sec  ',el4.6) 
if (y (2) . le.travp)go  to  303 

df ract- ( travp-disto) / (y ( 2 ) -disto) 
rravel-(y (1) -rmvelo) *dfract+rmvelo 
tmvel=(y (3) -tmvelo) *dfract+tmvelo 
write (3,318) rmvel , travel 

318  format (lx, 'muzzle  VELOCITY  m/s  ',el4.6,'  time  of  muzzle  velocity  s 

&ec  ' ,el4.6) 
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goto  319 

303  write(3,327)y(l) ,y(3) 

327  format (lx, 'velocity  of  projectile  m/s  #»el4.6,'  at  this  time  msec 
&' ,el4.6) 

319  efi=chwi*forcig/(gamai-l. ) 
efp=0 . 0 

do  315  i*l , nprop 

efp=efp+chwp(i) *forcp(i}/ (gamap(i)-l.O) 

315  continue 

tenerg=ef i+ef  p 
write (3,317) tenerg 

317  format (lx, 'total  initial  energy  available  J  =  ',el4.6) 
tenga  s=chwi *  f ore ig* tgas/ (garoai-1. ) /tempi 
do  135  i=l, nprop 

tengas=(frac(i)  *chwp(i)  *forcp(i)  *tgas/tempp(i)/(gamap(i)-l.)  )+teng 
&as 

write (3, 328) i,frac(i) 

328  format ('  FOR  PROPELLANT  ',12,'  MASSFRACT  BURNT  IS  ',el4.6) 

135  continue 

write (3 , 136) tengas 

136  format (lx, 'total  energy  remaining  in  gas  J=  ',el4.6) 
write (3,320) elpt 

320  format (lx, 'energy  loss  from  projectile  translation  J*=  ',el4.6) 
write (3 , 321) elpr 

321  format (lx, 'energy  loss  from  projectile  rotation  J=  ',el4.6) 
write ( 3 , 322 ) elgpm 

322  format (lx, 'energy  lost  to  gas  and  propellant  motion  J=  ',el4.6) 
write(3,323)elbr 

323  format (lx, 'energy  lost  to  bore  resistance  J=  ',el4.6) 
write(3,324)elrc 

324  format (lx, 'energy  lost  to  recoil  J=  ',el4.6) 
write(3,325)elht 

325  format (lx, 'energy  loss  from  heat  transfer  J-  ',el4.6) 
write (3 , 326) elar 

326  format (lx, 'energy  lost  to  air  resistance  j=  ',el4.6) 

c  call  gettim(ihro, imino, iseco, ihunso) 

c  time= ( ihro-ihr) *3600 . + ( imino-imin) *60 . + ( iseco-isec) + ( ihunso- ihuns) 

c  b/100. 

c  write (3,*) time 

stop 

20  write(*,140) 

140  format (lx, 'end  of  file  encounter') 
stop 

30  write(*,150) 

999  continue 

998  continue 

150  format (lx, 'read  or  write  error') 

stop 
end 

SUBROUTINE  PRF017 (P, PI, D, D1 , L, SURF, MASSF, X, NP, u) 

IMPLICIT  REAL*4(A-Z) 

C 

C  P=OUTEK  PERF  DIA 

C  P1=INNER  PERF  DIA 

C  D=OUTER  DIA 
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C  D1=DISTANCE  BETWEEN  PERF  CENTRES 
C  L-GRAIN  LENGTH 

C  NP*NUMBER  OF  PERFS 

C 

C  SURF-OUTPUT  SURFACE  AREA 

C  MASSF-OUTPUT  MASS  FRACTION  OF  PROPELLANT  BURNER 

C 

C  W*WEB  BETWEEN  OUTER  PERFS 

C  WO-OUTER  WEB 

C  Wl-WEB  BETWEEN  OUTER  AND  INNER  PERFS 

C  W4-MINIMUM  WEB 

C 

INTEGER  ITYM,NP 

DATA  PI ,SQRT3/3 . 14159 , 1 .732051/ ,  ITYtyO/ 

DATA  HAFPI,PIFOR,TWOPI/1.  570796, .785398,6.283185/ 
C 

IF(ITYM.GT.O)GO  TO  10 
P1SQ“P1*P1 
D1SQ=D1*D1 
PSQ=P*P 
DSQ=D*D 
D1SQ3«D1*SQRT3 
D2SQ3=D1SQ*SQRT3 
IF(NP.EQ.O)GO  TO  2000 
IF(NP.EQ.l)GO  TO  3000 
IF(NP.NE.7)GO  TO  60 
IF (PI .GT. (PfDl* (SQRT3-1) ) )  GO  TO  60 
IF(D.GE. Dl* (SQRT3+1. ) -P)GO  TO  130 
60  WRITE (6, 90) 

9  0  FORMAT  (IX, '  UNACCEPTABLE  GRANULATION ' ) 

STOP 

130  W=D1-P 

IF (W. LT. 0)GO  TO  60 
W0=(D-P-2.*Dl)/2. 

IF(W0.LT. 0. )GO  TO  60 
Wl=(2.*Dl-P-Pl)/2. 

IF(W1.LT. 0. )GO  TO  60 

X1=(P1SQ-PS<H4. *DlSQ-2 . *Pl*DlSQ3)/4 ./ (D1SQ3+P-P1) 
X2“  (4 . *DlSQ+D*D-2 . *D*D1SQ3-PSQ) /4 ./ (-D1SQ3+P+D) 
A=PI*L*  ( D+Pl+6 . *P) +HAFPI* (D  71SQ-6 . *PSQ) 

U=PI*I/4 . * (DSQ-P1SQ-6 . *PSQ) 

W4»AMIN1(W,W0,W1) 

10  MASSF-0. 

TWOX=X+X 

XSQ=X*X 

P1P2X=P1+TW0X 

PP2X=P4-TWOX 

DM2X=D-TWOX 

LM2  X=L-TWOX 

P12XSQ=P1P2X*P1P2X 

PP2XSQ=PP2X*PP2X 

DM2XSQ=DM2X*DM2X 

IF(NP.EQ. 0)GO  TO  2000 

IF(NP.EQ. l)GO  TO  3000 

IF(IM2X.GT.O)GO  TO  340 
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SURF=Q . 

V=0. 

GO  TO  850 

340  S0=PI*LM2X* (D+Pl+6 . 12 . *X) +HAFPI* (EM2X*DM2X 
1  -PlP2X*PlP2X-6 . *PP2X*PP2X) 

VO=PIFOR*LM2X* (DM2X*DM2X-P1P2X*P1P2X~6 . *PP2X*PP2X) 
IF(X.GT.W4/2.)GO  TO  360 
MASSF=-TWOX/L/ (DSQ-P1SQ-6 . *PSQ) 

MASSF=MASSF* (24 . *XSQ+ (24 . *P+4 . *Pl+4 . *D-12 . *L) *X+P1SQ 
1  +6.*PSQ-2.*L*D-2.*P1*L-12.*L*P-DSQ) 

SURF=S0 

RETURN 

360  IF (X.GT.W1/2 . )GO  TO  390 
F2=0. 

L2=0. 

A3=0. 

A4=0 . 

GO  TO  460 

390  Z=(2.*Dl+P+Pl+4.*X)/4. 

B3=( ( Pl-P) * ( Pl+P+4 . *X) +4 . *D1SQ) /4 . /D1/P1P2X 
A3  =ATAN ( SQRT (l.-B3*B3)/B3) 

B4= ( (P-Pl) ♦ (P+Pl+4 . *X) +4 . *D1SQ) /4 ./D1/PP2X 
A4=ATAN (SQRT ( 1 . -B4  *B4 ) /B4) 

F2=A3/4 . *P12XSQ+A4/4 . *PP2XSQ 
1  -SQRT (Z* (Z-Dl) * (2 . *Z-P-TWOX) * (2 . *Z-Pl-TWOX) ) 
L2=LM2X* (A4*PP2X+A3*P1P2X) 

460  IF(X.GT.W/2. )GO  TO  490 
F3=0 . 

L3=0 . 

A5=0. 

GO  TO  530 
490  B5=D1/PP2X 

A5=ATAN (SQRT ( 1 . -B5*B5) /B5) 

F3=(A5* PP2 XSQ-D1 *SQRT ( PP2XSQ-D1SQ) ) / 2 . 

13-2 .  *A5*IW2X*PP2X 
530  IF (X .GT. WO/2 . )GO  TO  560 
F1=0 . 

L1=0 . 

A1=0. 

A2=0. 

GO  TO  650 

560  Y=(2.*Dl+P+D)/4. 

Bl= ( (D+P) * (D-P-4 . *X) -4 . *D13Q)/4 . /D1/PP2X 
A1=ATAN(SQRT(1.-B1*B1)/B1) 

IF (A1 . GT. 0 . ) GO  TO  610 
A1=PI+A1 

610  B2=  (  (D+P)  *  (D-P-4  .  *X)  +4  .  *DlSQ)/4  ./D1/DM2X 
A2=ATAN (SQRT ( 1 . -B2  *B2 ) /B2 ) 

Fl=Al/4 . *PP2XSQ-A2/4 . *DM2XSQ+SQRT ( Y* (Y-Dl) 

1  * (2 . *Y-P-TWOX) * (2 . *Y-D+TWOX) ) 
L1=LM2X*(A1*PP2X+A2*EM2X; 

650  IF(X.GT.W/2. )GO  TO  690 

SURF- S 0+1 2 . * ( F1+F2+F3 ) -6 . * (L1+L2+L3) 

V=V0+6 . * (F1+F2+F3) *LM2X 


690  IF(X.LT.X1)G0  TO  730 

51- 0.0 
Vl-0.0 

GO  TO  760 

730  S1«3.*D2SQ3-PI*PP2XSQ-HAFPI*P12XSQ 
$  +6.*F3+12.*F2 

S1*S1+IM2X* (2 . * (PI-3 . *A5-3 . *A4) *PP2X+ (PI-6 .  *A3 ) 

$  *P1P2X) 

V1-I242X/2 .  *  (3 .  *D2SQ3-PI*PP2XSQ 
$  -HAFPI*P12XSQ+6.*F3+12.*F2) 

760  IF(X.LT.X2)  GO  TO  800 

52- 0.0 
V2-0.0 

GO  TO  830 

800  S2»HAFPI*DM2XSQ-3.*D2SQ3-TWOPI*PP2XSQ 
$  +12.*F1+6.*F3 

S2=S2+IM2X* ( (PI-6.  *A2) *DM2X+2 . * (TWOPI-3 .  *Al-3 .  *A5) 

$  *PP2X) 

V2=I242X/2 .  *  (HAFPI*EM2XSQ-3  .  *D2SQ3-1W0PI 
$  *PP2XSQ+12.*F1+6.*F3) 

830  SURF=S1+S2 
V=V1+V2 

850  MASSF-l.-V/U 
RETURN 
C 

C  ZERO  PERF  CALCULATIONS  START  HERE. 

C 

2000  if (d-2*x.le.0.0)  go  to  2001 
twox=x+x 
xsq=x*x 

MASSF=TWOX*  ( DSQ+  2 .  *L*D-4 .  *X*D-TWOX*L+4 .  *XSQ)  /  (  DSQ*L) 
U=dsq*l*pi/4. 

SURF=PI* (DSQ/2 . -4 . *D*X-TWOX*LfD*Lf  6 . *XSQ) 

RETURN 

2001  surf  *=0.0 
mass f= 1.0 
u=dsq*l*pi/4 . 
return 
C 

C  ONE  PERF  CALCULATIONS  START  HERE. 

C 

3000  if (d-p-4.*x.le.0.0)  goto  3001 
twox=x+x 

MASSF=TY70X*  ( DSQ+2 .  *L*D-4  .  *X*D-PSQ+  2  .  *P*L-4  .  *P*X) 

$  / (DSQ*L-P5Q*L) 
u=dsq*l*pi/4 . -psq*l*pi/4 . 

SURF=PI*  (DSQ/2  .  -4 .  *D*X-4  .  *X*P+D*L*P*L-PSQ/2 . ) 

RETURN 

3001  surf =0.0 
mass f= 1.0 

u=dsq*l*pi/4 . -psq*l*pi/4 . 

return 

END 
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APPENDIX  C 


USER’S  MANUAL  FOR  IBRGA 

IBRGA  relies  on  an  input  data  base  consisting  of  all  numerical  parameters  essential  for 
running  the  code.  All  values  are  in  metric  units.  Below  is  a  compilation  of  a  typical  IBRGAC 
data  base  showing  the  name  and  location  of  each  input  parameter.  The  names  for  the 
numerical  values  are  prefixed  with  an  alphabetical  designator  corresponding  to  the  position  at 
which  the  data  is  to  appear,  that  is,  from  left  to  right.  The  data  may  be  separated  by  blanks  or 
commas.  The  units  are  shown  to  the  right  of  each  input. 

ABCDEFGHIJK 
record  1 

A.  -  chamber  volume  (cm3) 

E.  -  groove  diameter  (cm) 

C.  -  land  diameter  (cm) 

D.  -  groove/land  ratio  (none) 

E.  -  twist  (turns/caliber) 

F.  -  projectile  travel  (cm) 

G.  -  gradient  switch  (  1  «  Lagrange,  2  =  chambrage  ) 

record  la  (  Read  only  if  gradient  switch  =  2  ) 

A.  -  number  of  points  to  describe  chamber  (/  <  10  ) 

B.  -  initial  distance  from  breech  (  must  be  0.0  cm  ) 

C.  -  diameter  at  0  (cm) 


Ith  distance  from  breech  (  position  of  base  of  projectile  (cm)) 
Ith  diameter  at  Ith  distance  (  used  to  calculate  bore  area  ).(cm) 

record  2 

A.  -  projectile  mass  (kg) 

B.  *  switch  to  calculate  energy  lost  to  nir  resistance 

C.  -  fraction  of  work  done  against  bore  to  heat  tube 

D.  -  gas  pressure  in  front  of  projectile  (Pa) 

record  3 

A.  -  number  of  barrel  resistance  points  (/  <  10) 

B.  -  bore  resistance  (MPa) 

C.  -  travel  (cm) 


Jth  bore  resistance  (MPa) 
Jth  travel  (cm) 
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record  4 

A.  -  mass  of  recoiling  parts  (kg) 

B.  -  number  of  recoil  point  pairs  (  2  ) 

C.  *  recoil  force  (N) 

D.  -  recoil  time  (s) 

E.  -  recoil  force  (N) 

F.  -  recoil  time  (s) 

record  5 

A.  -  free  convective  heat  transfer  coefficient  (  W/cm2  ) 

B.  -  chamber  wall  thickness  (cm) 

G  -  heat  capacity  of  steel  of  chamber  wall  (J/g-K) 

D.  -  initial  temperature  of  chamber  wall  (K) 

E.  -  heat  loss  coefficient 

F.  -  density  of  chamber  wall  steel  (g/cm3) 


record  6 

A.  -  impetus  of  igniter  propellant  (J/g) 

B.  -  covolume  of  igniter  (cm3/g) 

G  -  adiabatic  flame  temperature  of  igniter  propellant  (K) 

D.  -  initial  mass  of  igniter  (kg) 

E.  -  ratio  of  specific  heats  for  igniter 

record  7 

A.  -  number  of  propellants  (  K  <  10) 

B.  -  impetus  of  propellant  (J/g) 

G  -  adiabatic  temperature  of  propellant  (K) 

D.  -  covolume  of  propellant  (cm3/g) 

E.  -  initial  mass  of  propellant  (kg) 

F.  •  density  of  propellant  (g/cm3) 

G.  -  ratio  of  specific  heats  for  propellant 

H.  -  number  of  perforations  of  propellant  (  may  be  0,1,  or  7  only  ) 

I.  -  length  of  propellant  grain  (cm) 

J.  -  diameter  of  inner  perforations  in  propellant  grains  (cm) 

K.  -  diameter  of  outer  perforations  of  propellant  grains  (cm) 

L  -  outside  diameter  of  propellant  gram  (cm) 

M.  -  distance  between  perf  centers  (cm) 


(  Kth  propellent  ) 

A.  -  impetus  of  propellant  (J/g) 

B.  -  adiabatic  temperature  of  propellant  (K) 

G  -  covolume  of  propellant  (cm3/g) 

D.  -  initial  mass  of  propellant  (kg) 

E.  -  density  of  propellant  (gm/ern  ) 

F.  -  rat'O  of  specific  heats  for  propellant 

G.  -  number  of  perforations  of  propellant 

H.  -  length  of  propellant  grain  (cm) 

I.  -  diameter  of  inner  perforations  in  propellant  grains  (cm) 

J.  -  diameter  of  outer  perforations  of  propellant  grains  (cm) 

K.  -  outside  diameter  of  propellant  grain  (cm) 

L.  -  distance  between  perf  centers  (cm) 
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record  8 

A.  -  number  of  burning  rate  points  ( /  <  10  )  for  propellent  1 

B.  -  exponent 

C.  -  coefficient  (cm/s -MPa01  ) 

D.  -  pressure  (MPa) 


Jth  exponent 
Jth  coefficient 
Jth  pressure 


A.  -  number  of  burning  rate  points  (  L  <  10  )  for  propellent  N 

B.  -  exponent 

C.  -  coefficient  (cm/s -MPa***  ) 

D.  -  pressure  (MPa) 


Lth  exponent 
Lth  coefficient 
Lth  pressure 

record  9 

A.  -  time  increment  (ms) 

B.  -  print  increment  (ms) 

C  -  time  to  stop  calculation  (ms) 


APPENDIX  D 


Input  data  base  for  the  bo re -diameter  calculation. 


9832.2384  12.7  12.7  1.0  0.0  457.2  1 

9.796  0  0.0  0.0 

5  0.0  0.0  0.0  .6  0.0  1.3  0.0  300.  0.  457. 

I.e20  2  3.0e+4  0.0  8.0e+5  0.2 
0.001135  .01143  .46028  273.  1.  7.8612 
84.5535  .9755  294.  .004712  1.4 

1  1135.99  3141.  .9755  8.7  1.6605  1.23  7  3.175  .0508  .0508  1.07208908  .28072226 

1  1.0  .1105187  689.476 
0.005  .05  30. 
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APPENDIX  D 


Output  for  the  bo re -diameter  calculation. 

THE  INPUT  FILE  IS  iclopt 
using  Lagrange  pressure  gradient 
chamber  vclume  cm**3  0.983224E+04 

groove  diam  cm  0.127000E+Q2 
land  diam  cm  0.127000E+02 
groove/land  ratio  0.100000E+01 
twist  turns/ caliber  O.OOOOOOE+OO 
projectile  travel  cm  0.457200E+03 


projectile  mass  kg  0.979600E+01 
switch  to  calculate  energy  lost  to  air  resistance  J 
fraction  of  work  against  bore  used  to  heat  the  tube 
gas  pressure  Pa  O.OOOOOOE+OO 
number  barrel  resistance  points  5 
bore  resistance  MPa  -  travel  cm 


O.OOOOOOE+OO 
0 . 000000E+00 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


O.OOOOOOE+OO 
0.600000E+00 
0 . 130000E+01 
0.300000E+03 
0.457000E+03 


0 

O.OOOOOOE+OO 


mass  of  recoiling  parts  kg  0.100000E+21 
number  of  recoil  point  pairs  2 
recoil  force  N  recoil  time  sec 
0 . 300000E+05  0. OOOOOOE+OO 

0.800000E+06  0. 200000E+00 


free  convective  heat  transfer  coefficient  w/cm**2  K  0. 113500E-02 
chamber  wall  thickness  cm  0.114300E-01 

heat  capacity  of  steel  of  chamber  wall  J/g  K  0.460280E+00 
initial  temperature  of  chamber  wall  K  0.273000E+03 
heat  loss  coefficient  0.100000E+01 
density  of  chamber  wall  steel  g/cm**3  0.786120E+01 


impetus  of  igniter  propellant  J/g  0.845535E+02 
covolume  of  igniter  cm**3/g  0.975500E+00 

adiabatic  flame  temperature  of  igniter  propellant  K  0.294000E+03 
initial  mass  of  igniter  kg  0.471200E-02 
ratio  of  specific  heats  for  igniter  0.140000E+01 


for  propellant  number  1 
impetus  of  propellant  J/g  0.113599E+04 
adiabatic  temperature  of  propellant  K  0.314100E+04 
covolume  of  propellant  cm**3/g  0.975500E+00 

initial  mass  of  propellant  kq  0.870000E+01 
density  of  propellant  g/cm**3  0.166050E+01 

ratio  of  specific  heats  for  propellant  0.123000E+01 
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number  of  perforations  of  propellant  7 
length  of  propellant  grain  cm  0.317500E+01 

diameter  of  inner  perforation  in  propellant  grains  cm  0.508000E-01 
diameter  of  outerperf oration  of  propellant  grains  cm  0.508000E-01 
outside  diameter  of  propellant  grain  cm  0.107209E+01 
distance  between  perf  centers  cm  0.28O722E+C0 


number  of  burning  rate  points  1 
exponent  coefficient 

cm/sec-MPa**ai 

0. 100000E+Q1  0. 110519E+00 


pressure 

MPa 

0. 689476E+03 


time  increment  msec  0.500000E-02  print  increment  msec  0.500000E-01 
time  to  stop  calculation  msec  0.300000E+02 

area  bore  mA2  0.126677E-01  pressure  from  ign  Pa  0.868339E+05 
volume  of  unbumt  prop  mA3  0.523939E-02 

init  cham  vol-cov  ign  m  A3  0.982764E-02 

time  acc  vel  dis  mpress  pbase  pbrch 

0 . 5000E-05  0.8761E+02  0.435SE-03  0.1087E-08  0.8782E+05  0.6775E+05  0.9785E+05 
0 . 5000E-04  0. 9693E+02  0.4585E-02  0.1125E-06  0.9716E+05  0.7496E+05  0.1083E+06 

0 . 1050E-03  0. 1095E+03  0.1026E-01  0.5174E-06  0.1098E+06  0.8471E+05  0.1223E+06 

0 . 1550E-03  0.1223E+03  0.1605E-01  0.1172E-05  0.1226E+06  0.9456E+05  0.1366E+06 

0 • 2000E-03  0. 1349E+03  0.2183E-01  0.2022E-05  0.1352E+06  0.1043E+06  0.1506E+06 

0 . 2500E-03  0. 1502E+03  0.2895E-01  0.3289E-05  0.1505E+06  0.1161E+06  0.1677E+06 

0 . 3050E—03  0.1687E+03  0.3771E-01  0.5117E-05  0.1691E+06  0.1305E+06  0.1885E+06 

0 . 3550E-03  0. 1873E+03  0.4660E-01  0.7221E-05  0.1878E+06  0.1449E+06  0.2092E+06 

0 . 4050E-03  0.2077E+03  0.5647E-01  0.9793E-05  0.2081E+06  0.1606E+06  0.2319E+06 

0.4550E-03  0.2298E+03  0.6740E-01  0.1289E-04  0.2304E+06  0.1777E+06  0.2567E+06 

0 . 5050E-03  0. 2539E+03  0.7949E-01  0.1655E-04  0.2545E+06  0.1964E+06  0.2836E+06 

0 . 5550E— 03  0. 2801E+03  0.9283E-01  0.2086E-04  0.2807E+06  0.2166E+06  0.3128E+06 

0 . 6000E-03  0.3055E+03  0.1060E+00  0.2532E-04  0.306?E+Ofi  0.2762E+O6  0.3412E+06 

0 . 6500E-Q3  0.3358E+03  0.1220E+00  0.3102E-04  0.3366E+06  0.2597E+06  0.3751E+06 

0 . 7000E-03  0. 3686E+03  0.1396E+00  0.3755E-04  0.3695E+06  0.2850E+06  0.4117E+06 

0.7500E-03  0. 4039E+03  0.1589E+00  0.4501E-04  0.4048E+06  0.3123E+06  0.4511E+06 

0.8000E-03  0. 4417E+03  0.1800E+00  0.5347E-04  0.4428E+06  0.3416E+06  0.4934E+06 

0 . 8500E-03  0. 4824E+03  0.2031E+00  0.6305E-04  0.4835E+06  0.3730E+06  0.5388E+06 

0 . 9000E-03  0.5260E+03  0.2283E+00  0.7382E-04  0.5272E+06  0.4067E+06  0.5874E+06 

0 . 9500E-03  0. 5726E+03  0.2558E+00  0.8592E-04  0.5739E+06  0.4428E+06  0.6395E+06 

0 . 1000E-02  0. 6224E+03  0.2856E+00  0.9944E-04  C.6239E+06  0.4813E+06  0.6951E+06 

0 . 1050E-02  0. 6756E+03  0.3181E+00  0.1145E-03  0.6772E+06  0.5225E+06  0.7546E+06 

0 . 1100E-02  0. 7324E+03  0.3533E+00  0.1313E-03  0.7341E+06  0.5664E+06  0.8180E+06 

0.1150E-02  0. 7930E+03  0.3914E+00  0.1499E-03  0.7949E+06  0.6132E+06  0.8857E+06 

0.1200E-02  0. 8576E+03  0.4326E+00  0.1705E-03  0.8596E4-06  0.6632E+06  0.9579E+06 

0.1250E-02  0. 9265E+03  0.4772E+00  0.1932E-03  0.9287E-I-06  0.7165E+06  0.1035E+07 

0 . 1300E-02  0. 9999E+03  0.5254E+00  0.2183E-03  0.1002E+07  0.7733E+06  0.1117E+07 

0 . 1350E-02  0. 1078E+04  0.5773E+00  0.2458E-03  0.1081E+07  0.8338E+06  0.1204E+07 

0 . 1400E-02  0 . 1162E+04  0.6333E+00  0.2761E-03  0.1164E+07  0.8983E+06  0.1297E+07 

0.1450E-02  0. 1251E+04  0.6935E+00  0.3092E-03  0.1253E+07  0.9670E+06  0.1397E+07 

0 . 1500E-02  0 . 1345E+04  0.7584E+00  0.3455E-03  0.1348E+07  0.1040E+07  0.1503E+07 

0 . 1550E-02  0. 1446E+04  0.8282E+00  0.3851E-03  0.I450E+07  0.1119E+07  0.1615E+07 

0. 1600E-02  9. 1554E+04  0.9032E+00  0.4284E-03  0.1558E+07  0.1202E+07  0.1736E+07 

0 . 1650E-02  0 . 1669E-*-04  0.9837E+00  0.4755E-03  0,1673E-t-07  0.1291E+07  0.1864E+07 

0 . 1700E-02  0 . 1792E+04  0.1070E+01  0.5269E-03  0.1796E+07  0. 1386E+07  0.2001E+07 

0 . 1750E-02  0.1923E+04  0.1163E+01  0.5827E-03  0.1927E+07  0.1487E+07  0.2147E-K07 


J 


« 
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0. 1800E-Q2  0.2062E4-04 
0.1850E-02  0. 2211E+04 
0.1900E-02  0 . 2370E+04 
0. 1950E-02  0. 2540E+04 
0.2000E-02  0.2721E+04 
0 . 2C50E-02  0.2915E+04 
0 . 2100E-02  0 . 3122E+04 
0.2150E-02  0 . 3342E+04 
0 . 2200E-02  0.3578E+04 
0.2255E-02  0.3856E+04 
0.2305E-02  0.4127E+04 
0 . 2355E-02  0.4416E+04 
0.2405E-02  0. 4725E+04 
0.2455E-02  0.5055E+04 
0 . 2505E-02  0 . 5407E+04 
0 . 2555E-02  0.5783E+04 
0 . 2605E-02  0. 6185E+04 
0 • 2655E-02  0.6614E+04 
0 . 2705E-02  0. 7072E+04 
0.2755E-02  0.7561E+04 
0.2805E-02  0. 8083E+04 
0.2855E-02  0.8640E+04 
0 . 2905E-02  0. 9235E+04 
0.2955E-02  0.9869E+04 
0 . 3005E-02  0. 1055E+05 
0 . 3055E-02  0. 1127E+05 
0 . 3105E-02  0.1204E+05 
0 . 3155E-02  0. 1286E+05 
0 . 3205E-02  0. 1373E+05 
0 . 3255E-02  0.1466E+05 
0 . 3305E-02  0. 1566E+05 
0.3355E-02  0.1672E+05 
0 . 3405E-02  0 . 1784E+05 
0 . 3455E-02  0.1904E+05 
0 . 3505E-02  0. 2031E+05 
0 . 3555E-02  0.2167E+05 
0 . 3605E-02  0.2311E+05 
0 . 3655E-02  0.2464E+05 
0 . 3705E-02  0. 2627E+05 
0 . 3755E-02  0.2799E+05 
0 . 3805E-02  0. 2982E+05 
0 . 3855E-02  0.3177E+05 
0 . 3905E-02  0. 3382E+05 
0.3955E-02  0.3600E+05 
0 .  4005E-02  0. 3831E+05 
0 . 4055E-02  0. 4075E+05 
0 . 4 105E-02  0. 4332E+05 
0.4155E-02  0 . 4604E+05 
0 . 4205E-02  0. 4891E+05 
0.4255E-02  0.5193E+05 
0 . 4305E-02  0 . 5512E+05 
0 . 4355E-02  0.5846E+05 

r\  a  a  r\  cr  r*  a  r  nor  i 
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0 . 4 455E-02  0.6567E+05 
0 . 4 505E-02  0 . 6954E+05 


0. 1263E+01  0 . 6433E-G3 
0. 1369E+01  0.7091E-03 
0. 1484E+01  0.7804E-03 
0.1607E+01  0. 8576E-03 
0. 1738E+01  0.9412E-03 
0 . 1879E+01  0. 1032E-02 
0.2030E+01  0. 1129E-02 
0.2191E+01  0. 1235E-02 
0.2364E+01  0 . 1349E-02 
0.2569E+01  0 . 1484E-02 
0.2768E+01  0. 1617E-02 
0.2982E+01  0. 1761E-02 
0.3210E+01  0. 1916E-02 
0. 3454E+01  0. 2082E-02 
0. 3716E+01  0.2262E-02 
0.3996E+01  0.2454E-02 
0 . 4295E+01  0. 2662E-02 
0 . 4615E+01  0. 2884E-02 
0.4957E+01  0 . 3123E-02 
0.5322E401  0 . 3380E-02 
0.5713E+01  0.3656E-02 
0.6131E+01  0. 3952E-02 
0.6578E+01  0. 4270E-02 
0.7055E+01  0. 4610E-02 
0.7566E+01  0.4976E-02 
0.8111E+01  0. 5367E-02 
0.8693E+01  0.5787E-02 
0.9315E+01  0. 6237E-02 
0.9980E+01  0. 6720E-02 
0 . 1069E+02  0.7236E-02 
0. 1145E+02  0.7789E-02 
0. 1226E+02  0. 8382E-02 
0. 1312E+02  0.9016E-02 
0. 1404E+02  0.9695E-02 
0 . 1503E+02  0.1042E-01 
0. 1607E+02  0. 1120E-01 
0. 1719E+02  0.1203E-01 
0 . 1839E+02  0. 1292E-01 
0 . 1966E+02  0. 1387E-01 
0. 21Q2E+02  0.1489E-01 
0 . 2246E+02  0.1597E-01 
0 . 2400E+02  0. 1713E-01 
0 . 2564E+02  0.1837E-01 
0 . 2738E+02  0.1970E-01 
0 . 2S24E+02  0. 2111E-01 
0 . 3122E+02  0 . 2263E-01 
0 . 3332E+02  0.2424E-01 
0 . 3555E+02  0.2596E-01 
0 . 3792E+02  0.2780E-01 
0.4044E+02  0. 2975E-01 
0 . 4312E+02  0.3184E-01 
0 . 4596E+02  0.3407E-01 

i\  Ane\nx*\r\*\  a  a  c  a  a  t?  __  a 
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0.5216E+02  0. 3397E-01 
0 . 5554E+02  0.4166E-01 


0.2067E+07  0. 1595E+07 
0. 2216E+07  0.1710E+07 
0.2376E+07  0.1833E+07 
0.2546E+07  0.1964E+07 
0.2728E+07  0.2104E+07 
0 . 2922E+07  0.2254E+07 
0 . 3129E+07  0. 2414E+07 
0 . 3350E+07  0. 2585E+07 
0. 3587E+07  0.2767E+07 
0.3865E+07  0.2982E+07 
0.4137E+07  0. 3191E+07 
0.4426E+07  0.3415E+07 
0.4736E+07  0.3654E+07 
0 . 5067E+07  0.3909E+07 
0 . 5420E+07  0. 4181E+07 
0 . 5797E+07  0.4472E+07 
0 . 6200E+07  0. 4783E+07 
0 . 6630E+07  0. 5115E+07 
0.7089E+07  0. 5469E+07 
0.7579E+07  0.5847E+07 
0.8102E+07  0. 6251E+07 
0.8661E+07  0. 6682E+07 
0.9257E+07  0.7141E+07 
0.9893E+07  0.7632E+07 
0 . 1057E+08  0.8155E+07 
0 . 1129E+08  0. 8714E+07 
0 . 1207E+08  0. 9309E+07 
0 . 1289E+08  0. 9943E+07 
0 . 1377E+08  0. 1062E+08 
0 .  l**70E+08  0 . 1134E+08 
0 . 1569E+08  0.1211E+0R 
0.1675E+08  0. 1293E+08 
0.1788E+08  0. 1380E+08 
0 . 1908E+08  0. 1472E+08 
0 . 2036E+08  0.X571E+08 
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0.2138E4-09 
0.2093E4-09 
0.2049E4-09 
0.2006E4-09 
0.19C4E4-09 
0 . 1923E+09 
0.1882E4-09 
0.1843E4-09 
0 . 1805E+09 
0.1768E409 
0.1732E+09 
0 : 1697EX09 
0.1663E+09 
0. 1630E+09 
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0.1005E-01  0. 1430E+06  0.J-317E+04  0.37Q7E+01  0.1433E+09  0.1106E+09  0.1597E+09 

0.1010E-01  0. 1402E+06  0.1324E+04  0.3773E+01  0.1405E+09  0.1084E+09  0.1566E+09 

0.1015E-01  0. 1374E+06  0.1331E-f04  0.3840E+01  0.1377E+09  0.1063E+09  0.1535E+09 

0. 1020E-01  0. 1347E+06  0.1338E+04  0.3906E+01  0.1351E+09  0.1042E+09  0.1505E+09 

0.1025E-01  0. 1321E+06  0.1345E+04  0.3974E+01  0.1325E+09  0.1022E+09  0.1476E+09 

0.1030E-01  0.1296E+06  0.1351E+04  0.4041E+01  0.1299E+09  0.1002E+09  0.1448E+09 

0. 1035E-01  0. 1271E+06  0.1358E+04  0.4109E+01  0.1274E+09  0.9832E+08  0.1420E+09 

0 . 1040E-01  0. 1247E+06  0.1364E+04  0.4177E+01  0.1250E+09  0.9646E+08  0.1393E+09 

0 . 1045E-01  0. 1224E+06  0.1370E+04  0.4245E+01  0.1227E+09  0.9466E+08  0.1367E+09 

0.1050E-01  0. 1201E+06  0.1376E+04  0.4314E+01  0.1204E+09  0.9289E+08  0.1342E+09 

0 . 1055E-01  0. 1179E+06  0.1382E+04  0.4383E+01  0.1182E+09  0.9118E+08  0.1317E+09 

0.1060E-01  0. 1158E+06  0.1388E+04  0.4452E+01  0.1160E+09  0.8951E+08  0.1293E+09 

0. 1065E-01  0. 1136E+06  0.1394E+04  0.4521E+01  0.1139E+09  0.8788E+08  0.1269E+09 

del tat  t  0.106950E-01  intg  t  0.106948E-01 
PMAXMEAN  Pa  0.310784E+09  time  at  PMAXMEAN  sec  0.707511E-02 

PMAXBASE  Pa  0.239766E+09  time  at  PMAXBASE  sec  0.707511E-02 

PMAXBREECH  Pa  0.346294E+09  time  at  PMAXBREECH  sec  0.707511E-02 
muzzle  VELOCITY  nv's  0.139773E+04  time  of  muzzle  velocity  sec  0.106910E-Q1 
total  initial  energy  available  J  =  0.429710E+08 

FOR  PROPELLANT  1  MASSFRACT  BURNT  IS  0.996610E+00 
total  energy  remaining  in  gas  J=  0. 28957 5E+08 

energy  loss  from  projectile  translation  0.957471E+07 

energy  loss  from  projectile  rotation  J=  O.OOOOOOE+OO 
energy  lost.  ;o  gas  and  propellant  motion  J=  0.283602E+07 
energy  ] ost  to  bore  resistance  J=  O.OOOOOOE+OO 
energy  io.'t  to  recoil  J=*  0.187099E-11 

energy  loss  from  heat  transfer  J=  0.145715E+07 
energy  lost  to  air  resistance  J=  O.OOOOOOE+OO 
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APPENDIX  E 


Input  data  base  for  the  chambrage  calculation. 

9832.2384  12.7  12.7  1.0  0.0  457.2  2 

3  0.  15.39110755  46.4820  15.39110755  54.1020  12.70 
9.796  0.  0.0  0.0 

5  0.0  0.0  0.0  .6  0.0  1.3  0.0  300.  0.  457. 

I.e20  2  3 .0e+4  0.0  8.0e+5  0.2 
.001135  .01143  .46028  273.  1.  7.8612 
84.5535  .9755  294.  .004712  1.4 

1  1135.99  3141.  .9755  8.85  1.6605  1.23  7  3.175  .0508  .0508  1.06290007  .27842504 
1  1.0  .1105187  689.476 
.005  .05  30. 
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APPENDIX  E 


Output  for  the  chambrage  chamber  calculation. 

THE  INPUT  FILE  IS  iclch2 

Using  chambrage  pressure  gradient 


chamber  distance  cm  chamber  diameter  cm 


0. OOOOOOE+OO  0. 153911E+02 

0.464820E+02  0.153911E+02 

0.541020E+02  0.127000E+02 

chamber  volume  cm**3  0.983215E+04 

groove  diam  cm  0.127000E+02 
land  diam  cm  0.127000E+02 
groove/land  ratio  0.100000E+01 
twist  tums/caliber  0. 000000E+00 
projectile  travel  cm  0.457200E+03 


projectile  mass  kg  0.979600E+01 
switch  to  calculate  energy  lost  to 
fraction  of  work  against  bore  used 
gas  pressure  Pa  0. 000000E+00 
number  barrel  resistance  points  5 
bore  resistance  MPa  -  travel  cm 


0. 000000E+00 
0. OOOOOOE+OO 
0. 000000E+00 
0. 000000E+00 
0. OOOOOOE+OO 


0. OOOOOOE+OO 
0. 600000E+00 
0. 130000E+01 
0 . 300000E+03 
0 . 457000E+03 


air  resistance  J  0 

to  heat  the  tube  0.000000E+00 


mass  of  recoiling  parts  kg  0.100000E+21 
number  of  recoil  point  pairs  2 
recoil  force  N  recoil  time  sec 
0 . 300000E+05  0. OOOOOOE+OO 

0 . 800000E+06  0. 200000E+00 


free  convective  heat  transfer  coefficient  w/cm**2  K  0.113500E-02 
chamber  wall  thickness  cm  0.114300E-01 

heat  capacity  of  steel  of  chamber  wall  J/g  K  0.460280E+00 
initial  temperature  of  chamber  wall  K  0.273000E+03 
heat  loss  coefficient  0.100000E+01 
density  of  chamber  wall  steel  g/cm**3  0.786120E+01 


impetus  of  igniter  propellant  J/g  0.845535E+02 
covolume  of  igniter  cm**3/g  0.975500E+00 

adiabatic  flame  temperature  of  igniter  propellant  K  0.294000E+03 
initial  mass  of  igniter  kg  0.471200E-02 
ratio  of  specific  heats  for  igniter  0.140000E+01 


for  propellant  number  1 
impetus  of  propellant  J/g  0.113599E+04 
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adiabatic  temperature  of  propellant  K  0.314100E+04 
covolume  of  propellant  cm**3/g  0.975500E+00 

initial  mass  of  propellant  kg  0.885000E+01 
density  of  propellant  g/cm**3  0.166050E+01 

ratio  of  specific  heats  for  propellant  0.123000E+01 
number  of  perforations  of  propellant  7 
length  of  propellant  grain  cm  0.317500E+01 

diameter  of  inner  perforation  in  propellant  grains  cm  0.508000E-01 
diameter  of  outerperforat.ion  of  propellant  grains  cm  0. 508000E-01 
outside  diameter  of  propellant  grain  cm  0.106290E+01 
distance  between  perf  centers  cm  0.278425E+00 


number  of  burning  rate  points  1 

exponent  coefficient  pressure 

cm/sec-MPa**ai  MPa 

0 . 100000E+01  0. 110519E+00  0.689476E+03 

time  increment  msec  0.500000F-02  print  increment  msec  0.500000E-01 
time  to  stop  calculation  msec  0.300000E+02 

area  bore  mA2  0.126677E-01  pressure  from  ign  Pa  0.885796E+05 
volume  of  unbumt  prop  mA3  0.532972E-02 

init  cham  vol-cov  ign  m  A3  0.982755E-02 

time  acc  vel  dis  mpress  phase  pbrch 

0 . 5000E-05  0.9985E+02  0.4141E-03  0.8272E-09  0.8963E+05  0.7721E+05  0.9509E+05 
0 . 5000E-04  0. 1110E+03  0.5154E-02  0.1242E-06  0.9964E+05  0.8584E+05  0.1057E+06 

0 . 1050E-03  0. 1262E+03  0.1167E-01  C.5830E-06  0.1133E+06  0.9757E+05  0.1202E+06 

0 . 1550E-03  0 . 1416E+03  0.1835E-01  0.1330E-05  0.1271E+06  0.1095E+06  0.1348E+06 

0 . 2000E-03  0. 1568E+03  0.2506E-01  0.2305E-05  0.1408E+06  0.1213E+06  0.1493E+C6 

0. 2500E-03  0 . 1755E+03  0.3336E-01  0.3761E-05  0.1575E+06  0.1357E+06  0.1671E+06 

0.3050E-03  0 . 1982E+03  0.4363E-01  0.5873E-05  0.1779E+06  0.1533E+06  0.1888E+06 

0.3550E-03  0 . 2210E+03  0.5410E-01  0.8311E-05  0.1984E+06  0.1709E+06  0.2105E+06 

0 . 4050E-03  0 . 2461E+03  0.6577E-01  0.1130E-04  0.2209E+06  0.1903E+06  0.2343E+06 

0. 4550E-03  0 . 2735E+03  0.7875E-01  0.1491E-04  0.2455E+06  0.2115E+06  0.2604E+06 

0. 5050E-03  0 . 3033E+03  0.9315E-01  0.1920E-04  0.2723E+06  0.2346E+06  0.2889E+06 

0 - 5550E-03  0 • 3359E+03  0.1091E+00  0.2425E-04  0.3015E+06  0.2597E+06  0.3199E+06 

0.6000E-03  0. 3675E+03  0.1249E+00  0.2951E-04  0.3299E+06  0.2842E+06  0.3500E+06 

0 . 6500E-03  0. 4054E+03  0.1443E+00  0.3623E-04  0.3640E+06  0.3135E+06  0.3861E+06 

0 . 7000E-03  0 . 4465E+03  0.1655E+00  0.4397E-04  0.4008E+06  0.3452E+06  0.4252E+06 

0 . 7500E-03  0 . 4907E+03  0.1889E+00  0.5282E-04  0.4405E+06  0.3795E+06  0.4673E+06 

0 . &000E-G3  0. 5384E+03  0.2147E+00  0.6290E-04  0.4833E+06  0.4163E+06  0.5127E+06 

0 . 8500E-03  0. 5896E+03  0.2428E+00  0.7433E-04  0.5293E+06  0.4559E+06  0.5615E+06 

0 • 9000E-03  0 . 6446E+03  0.2737E+00  0.8723E-04  0.5787E+06  0.4985E+06  0.6140E+06 

0 . 9500E-03  0 . 7036E+03  0.3074E+00  0.1017E-03  0.6317E+06  0.5441E+06  0.6702E+06 

0 . 1000E-02  0.76S9E403  0.3441E+-00  0.1180E-03  0.6885E+06  0.5930E+06  0.7304E+06 

0 . 1050E-02  0 . 8346E+03  0.3841E+00  0.1362E-03  0.7493E+06  0.6454E+06  0.7949E+06 

0 - 1100E-02  0 . 9070E403  0.4277E+00  0.1565E-03  0.8143E+06  0.7014E+06  0.8639E+06 

0 . 115QE-02  0 . 9845E+03  0.4749E+00  0.1790E-03  0.8839E+06  0.7613E+06  0.9377E+06 

0.1200E-02  0 . 1067E+04  0.5262E+00  0.2041E-03  0.9583E+06  0.8254E+06  0.1017E+07 

0 . 1250E-02  0 . 1156E+04  0.5818E+00  0.2317E-03  0.1038E407  0.3939E+06  0.1101E407 

0 » 1300E-02  0. 1251E+04  0-6419E+00  0.2623E-03  0.1123E+07  0.9671E+06  0.1191E+07 

0 . 1350E-02  0.1352E+04  0.7069E+00  0.2960E-03  0.1214E+07  0-1045E+07  0.1288E+07 

0 . 1400E-02  0 . 14 60E+04  0.7772E+00  0.3331E-03  0.1311E+U/  U.112yE+U/  U.i3yiE+0/ 

0 . 1450E-02  0 . 1576E+04  0.8531E+00  0.3738E-03  0.1415E+07  0.1219E+07  0.1501E+07 

0 . 1500E-02  0.1700E+04  0.9349E400  0.4185E-03  0,1526E+07  0.1314E+-07  0.1619E+07 
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0.1550E-02 
0. 1600E-02 
0.1650E-02 
0.1700E-02 
0.1750E-02 
0.1800E-02 
0 . 1850E-02 
0.1900E-02 
0 , 1950E-02 
0.2000E-02 
0.2050E-02 
0.2100E-02 
0.2150E-02 
0.2200E-02 
0 . 2255E-02 
0.2305E-02 
0 . 2355E-02 
0 . 2405E-02 
0 . 2455E-02 
0 . 2505E-02 
0 . 2555E-02 
0. 2605E-02 
0.2655E-02 
0 . 2705E-02 
0.2755E-02 
0. 2805E-02 
0. 2855E-02 
0. 2905E-02 
0. 2955E-02 
0. 3005E-02 
0. 3055E-02 
0. 3105E-02 
0.3155E-02 
0 . 3205E-02 
0. 3255E-02 
0. 3305E-02 
0. 3355E-02 
0.3405E-02 
0.3455E-02 
0.3505E-02 
0.3555E-02 
0 . 3505E-02 
0 . 3655E-02 
0 . 3705E-02 
0 . 3755E-02 
0 . 3805E-02 
0 . 3855E-02 
0 . 39Q5E-02 
0 . 3955E-02 
0 . 4005E-02 
0.4055E-02 
0 . 4105E-Q2 
0.4155E-02 
0 . 4205E-02 
0.4255E-02 


0. 1832E+04 
0. 1974E+04 
0.2125E+04 
0.2287E+04 
0. 2461E+04 
C.2647E+04 
0. 2846E+04 
0. 3059E+04 
0. 3287E+04 
0. 3531E+04 
0. 3793E+04 
0.4074E+04 
0.4374E+04 
0.4696E+04 
0. 5076E+04 
0. 5447E+04 
0.5845E+04 
0. 6272E+04 
0. 6728E+04 
0. 7217E+04 
0.7740E+04 
0.8300E+04 
0 . 8899E+04 
0.9540E+04 
0.1023E+05 
0. 1096E+05 
0.U74E+05 
0.1258E+05 
0. 1348E+05 
0. 1444E+05 
0. 1546E+05 
0. 1655E+05 
0 . 1772E+05 
0. 1896E+05 
0.2029E+05 
0 . 2170E+05 
0. 2321E+05 
0 . 2481E+05 
0 . 2652E+05 
0 . 2834E+05 
0 , 3027E+05 
0 . 3231E+05 
0 . 3449E+05 
0 . 3680E+05 
0 . 3924E+05 
0 .  4183E+05 
0 . 4457E+05 
0 . 4746E+05 
0 . 5051E+05 
0 . 5373E+05 
0. 5712E+05 
0. 6068E+Q5 
0.6443E+05 
0.6836E+05 
0.7247E+05 


0.1023E+01 
0 . 1118E+01 
0. 1221E+01 
0. 1331E+01 
0.1450E+01 
0 . 1577E+01 
0.1715E+01 
0.1862E+01 
0.2021E+01 
0.219 1E+01 
0.2374E+01 
0.2571E+01 
0.2782E+01 
0.3008E+01 
0.3277E+01 
0 . 3540E+01 
0 . 3822E+01 
0.4125E+01 
0 . 4450E+01 
0 . 4798E+01 
0 . 5172E+01 
0 . 5573E+01 
0.6003E+01 
0 . 6463E+01 
0 . 6957E+01 
0 . 7487E+01 
0 . 8054E+01 
0.8662E4-01 
0.9314E+01 
0. 1001E+02 
0. 1076E+02 
0. 1156E+02 
0. 1242E+02 
0.1333E+02 
0. 1431E+02 
0.1536E+02 
0. 1648E+02 
0. 1768E+02 
0 . 1897E+02 
0. 2034E+02 
0. 2180E+02 
0.2337E+02 
0 . 2504E+02 
0 . 2682E+02 
0 . 2872E+02 
0 . 3074E+02 
0 . 3290E+02 
0 . 3520E+02 
0.3765E+02 
0.4026E+02 
0.4303E+02 
0,4597F.-f02 
0.4910E-I-02 
0 . 5242E+02 
0 . 5594E : 02 


0.4674E-03 
0.5209E-03 
0.5794E-03 
0. 6431E-03 
0. 7126E-03 
0.7882E-03 
0.8705E-03 
0.9598E-03 
0.1057E-02 
0. 1162E-02 
0. 1276E-02 
0. 1400E-02 
0. 1533E-02 
0. 1678E-02 
0.1851E-02 
0.2021E-02 
0. 2205E-02 
0.2404E-02 
0. 2618E-02 
0. 2849E-02 
0.3098E-02 
0. 3367E-02 
0.3656E-02 
0. 3968E-02 
0.4303E-02 
0.4664E-02 
0.5052E-02 
0. 5470E-02 
0.5919E-02 
0. 6402E-02 
0.6921E-02 
0.7479E-02 
0.8078E-02 
0. 8721E-02 
0.9412E-02 
0. 1015E-01 
0. 1095E-01 
0. 1180E-01 
0. 1272E-01 
0 . 1370E-01 
0. 1476E-01 
0 . 1588E-01 
0.1709E-01 
0. 1839E-01 
0 . 1978E-01 
0 . 2126E-01 
0. 2285E-01 
0. 2456E-01 
0. 2638E-01 
0 . 2832E-01 
0. 3041E-01 
0.3263E-01 
0.  3501E-01 
0. 3754E-01 
0. 4025E-01 


0. 1645E+07 
0.1772E+07 
0.1909E+07 
0.2054E+07 
0.2210E+07 
0.2377E+07 
0 . 2556E+07 
0.2748E+07 
0.2953E+07 
0.3173E+07 
0.3408E+07 
0.3661E+07 
0.3931E+07 
0 . 4220E+07 
0.4563E+07 
0.4897E+07 
0.5256E+07 
0.5640E+07 
0 . 6051E+07 
0 . 6492E+07 
0 . 6963E+07 
0.7469E+07 
0 . 8010E+07 
0 . 8589E+07 
0 . 9208E+07 
0.9871E+07 
0 . 1058E+08 
0.1134E408 
0 . 1215E+08 
0 . 1302E+08 
0. 1395E+08 
0. 1494E+08 
0. 1600E+08 
0. 1713E+08 
0. 1833E+08 
0 . 1962E+08 
0 . 2099E+08 
0 . 2246E+08 
0.2402E+08 
0. 2568E+08 
0 . 2744E+03 
0.2932E+08 
0 . 3132E+08 
0.3344E408 
0 . 3569E+08 
0.3808E+08 
0 . 4061E+08 
0.4329E+08 
0 . 4613E+08 
0 . 4912E+08 
0.5228E+08 
0 . 5562E+08 
0 . 5913E+08 
0 . 6282E+08 
0 . 6670E+08 


0. 1417E+07 
0. 1526E+07 
0. 1643E+07 
0. 1769E+07 
0. 1903E+07 
0.2047E+07 
0. 2201E+07 
0. 2365E+07 
0.2542E+07 
0.2731E+07 
0. 2933E+07 
0. 3150E+07 
0. 3382E+07 
0. 3631E+07 
0. 3925E+07 
0. 4213E+07 
0. 4520E+07 
0. 4850E+07 
0. 5203E+07 
0. 5581E+07 
0. 5985E+07 
0. 6418E+07 
0.6882E+07 
0. 7377E+07 
0. 7908E+07 
0. 8475E+07 
0.9082E+07 
0.9731E+07 
0. 1042E+08 
0.1116E+08 
0. 1196E+08 
0 . 1280E+08 
0. 1370E+08 
0. 1466E+08 
0 . 1569E+08 
0 . 1678E+08 
0.1795E408 
0. 1919E+08 
0. 2051E+08 
0 . 2191E+08 
0. 2340E+08 
0. 2499E+08 
0 . 2667E+08 
0 . 2845E+08 
0. 3035E+08 
0 . 3235E+08 
0 . 3446E+08 
0 . 3670E+08 
0.3906E+08 
0. 4155E+08 
0 . 4417E+08 
0. 4693E+08 
0 . 4982E+08 
0.5286E+08 
0. 5604E+08 


0. 1745E+07 
0. 1880E+07 
0.2025E-K)7 
0.2179E+07 
0 „ 2345E+07 
0 . 2522E+07 
0.2712E+07 
0.2915E+07 
0.3133E+07 
0.3366E+07 
0.3616E+07 
0 . 3883E+07 
0.4170E+07 
0.4477E+07 
0.4840E+07 
0.5195E+07 
0.557 6E+07 
0 . 5983E+07 
0 . 6419E+07 
0 . 6887E+07 
0.7387E+07 
0.7923E+07 
0 . 8497E+07 
0.9111E+07 
0.9769E+07 
0 . 1047E+08 
0. 1122E+08 
0 . 1203E+08 
0 . 1289E+08 
0 . 1381E+08 
0 . 1480E+08 
0 . 1585E+08 
0 . 1697E+08 
0. 1817E+08 
0. 1945E+08 
0. 2081E+08 
0. 2227E+08 
0. 2382E+08 
0 . 2548E+08 
0 . 2724E+08 
0. 2911E+08 
0. 3111E+08 
0 . 3323E+08 
0. 3548E+08 
0. 3787E+08 
0 . 4040E+08 
0. 4309E+08 
0. 4593E+08 
0 . 4894E+08 
0 . 5212E+08 
0.5548E+08 
0 . 5902E+08 
0 . 6275E+08 
0 . 6668E+08 
0.7080E+08 
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0.4305E-02 
0.4355E-02 
0.4405E-02 
0.4455E-02 
0.4505E-02 
0.4555E-02 
0.4605E-02 
0.4655E-02 
0.4705E-02 
0.4755E-02 
0.4805E-02 
0 . 4855E-02 
0. 4905E-02 
0 . 4955E-02 
0 . 5005E-02 
0 . 5050E-02 
0 . 5100E-02 
0.5150E-02 
0 . 5200E-02 
0 . 5250E-02 
0 . 5300E-02 
0.5350E-02 
0 . 5400E-02 
0 . 5450E-02 
0 . 5500E-02 
0 . 5550E-02 
0 . 5600E-02 
0.5650E-02 
0 . 5700E-02 
0 . 5750E-02 
0 . 5800E-02 
0 . 5850E-02 
0 . 5900E-02 
0 . 5950E-02 
0 . 6000E-02 
0 . 6050E-02 
0 . 6100E-02 
0 . 6150E-02 
0 . 6200E-02 
0 . 6250E-02 
0 . 6300E-02 
0 . 6350E-02 
0 . 6400E-02 
0 . 6450E-02 
0 . 6500E-02 
0.6550E-02 
0 . 6600E-02 
0 . 6650E-02 
0 . 6700E-02 
0 . 6750E-02 
0 . 6800E-02 
0.6850E-02 
0 . 6900E-02 
0 . 6950E-02 
0 . 7000E-02 


0.7678E+05 
0.8127E+05 
0. 8596E+05 
0. 9083E+05 
0.9589E+05 
0. 1011E+06 
0. 1066E+06 
0. 1122E+06 
0 . 1180E+06 
0. 1239E+06 
0. 1300E+06 
0.1362E+06 
0.1426E+06 
0.1490E+06 
0. 1556E+06 
0 . 1616E+06 
0 . 1683E+06 
0. 1750E+06 
0. 1818E+06 
0. 1886E+06 
0. 1953E+06 
0. 2020E+06 
0.2087E+06 
0. 2153E+06 
0 . 2218E+06 
0. 2281E+06 
0.2343E+06 
0. 2404E+06 
0 . 2462E+06 
0.2519E+06 
0 . 2573E+06 
0. 2626E+06 
0. 2675E+06 
0. 2723E+06 
0.2767E+06 
0. 2809E+06 
0. 2848E+06 
0. 2885E+06 
0. 2918E+06 
0. 2949E+06 
0. 2977E+C6 
0. 3002E+06 
0. 3024E+06 
0. 3044E+06 
0. 3061E+06 
0 . 3075E+06 
0.3086E+06 
0. 3095E+06 
0. 3102E+06 
0 . 3106E+06 
0.3108E+06 
0. 3108E+06 
0. 3106E+06 
0.3102E+06 
0 . 3096E+06 


0. 5967E+02 
0. 6362E+02 
0.6780E+02 
0.7222E+02 
0.7689E+02 
0.8181E+02 
0.8700E+02 
0.9247E+02 
0.9822E+02 
0.1043E+03 
0.1106E+03 
0.1173E+03 
0. 1242E+03 
0.1315E+03 
0 . 1391E+03 
0. 1463E+03 
0 . 1545E+03 
0.1631E+03 
0 . 1720E+03 
0.1813E+03 
0 . 1909E+03 
0.2008E+03 
0.2111E+03 
0.2217E+03 
0.2326E+03 
0 . 2439E+03 
0.2554E+03 
0.2673E+03 
0 . 2794E+03 
0 . 2919E+03 
0. 3046E+03 
0 . 3176E+03 
0. 3309E+03 
0. 3444E+03 
0 . 3581E+03 
0. 3720E+03 
0 . 3862E+03 
0.4005E+03 
0 . 4150E+03 
0.4297E+03 
0 . 4445E+03 
0.4595E+03 
0. 4745E+03 
0.4897E+03 
0 . 5050E+03 
0. 5203E+03 
0.5357E+03 
0 . 5512E+03 
0 . 5667E+03 
0 . 5822E+03 
0 . 5977E+03 
0 . 6133E+03 
0. 6288E+03 
0 . 6443E+03 
0. 6598E+03 


0. 4314E-01 
0. 4622E-01 
0. 4951E-01 
0.5301E-01 
0. 5673E-01 
0.6070E-01 
0. 6492E-01 
0. 6940E-01 
0. 7417E-01 
0. 7923E-01 
0.8460E-01 
0. 9030E-Q1 
0.9633E-01 
0. 1027E+00 
0. 1095E+00 
0. 1159E+00 
0.1234E+00 
0. 1314E+00 
0.1397E+00 
0. 1486E+00 
0. 1579E+00 
0. 1677E+00 
0. 1780E+00 
0. 1888E+00 
0. 2001E+00 
0. 2121E+00 
0.2245E+00 
0. 2376E+00 
0. 2513E+00 
0. 2655E+00 
0. 2805E+00 
0. 2960E+00 
0. 3122E+00 
0. 3291E+00 
0. 3467E+00 
0. 3649E+00 
0. 3839E+00 
0 . 4035E+00 
0 . 4239E+00 
0 . 4451E+00 
0 . 4669E+00 
0.4895E+00 
0. 5129E+00 
0 . 5370E+00 
0. 5618E+00 
0 . 5875E+00 
0. 6139E+00 
0. 6410E+00 
0.6690E+00 
0 . 6977E+00 
0 . 7272E+00 
0 . 7575E+00 
0 . 7885E+00 
0 . 8204E+00 
0. 8530E+00 


0.7077E+08 
0 , 7504E+08 
0.7950E+08 
0.8415E+08 
0.8901E+08 
0.9406E+08 
0.9930E+08 
0. 1047E+09 
0. 1104E+09 
0. 1162E+09 
0.1221E+09 
0. 1233E+09 
0. 1346E+09 
0. 1410E+09 
0. 1476E+09 
0. 1536E+09 
0. 1604E+09 
0. 1672E+09 
0. 1741E+09 
0. 1811E+09 
0. 1880E+09 
0. 1950E+09 
0. 2019E+09 
0. 2088E+09 
0. 2156E+09 
0. 2223E+09 
0. 2289E+09 
0. 2353E+09 
0. 2416E+09 
0.2477E+09 
0. 2535E+09 
0.2592E+09 
0. 2646E+09 
0. 2698E+09 
0.2747E+09 
0.2793E+09 
0.2836E+09 
0.2877E+09 
0 . 2914E+09 
0.2949E+09 
0 . 2980E+09 
0.3008E+09 
0. 3034E+09 
0 . 3056E+09 
0.3076E+09 
0.3092E+C9 
0 . 3106E+09 
0 . 3117E+09 
0.3125E+09 
0 . 313 1E+09 
0 . 3 135E+09 
0.3136E+09 
0 . 3135E+09 
0.3132E+U9 
0 . 3 127E+09 


0.5937E+08 
0.6285E+08 
0.6647E+08 
0.7024E+08 
0.7415E+08 
0. 7821E+08 
0.8241E+08 
0.8675E+08 
0.9121E+08 
0. 9581E+08 
0.1005E+09 
0. 1053E+09 
0.1102E+09 
0, 1152E+09 
0. 1203E+09 
0.1249E+09 
0.1301E+09 
0. 1353E-r09 
0. 1406E+09 
0.1458E+09 
0. 1510E+09 
0. 1562E+09 
0. 1614E+09 
0. 1665E+09 
0. 1715E+09 
0. 1764E+09 
0.1812E+09 
0.1859E+09 
0. 1904Er09 
0. 1948E+09 
0. 1990E+09 
0. 2030E+09 
0.  2O69E+0S 
0. 2105E+09 
0.  2140E+09 
0. 2172E+09 
0. 2203E+09 
0.2231E+09 
0. 2257E+09 
0. 2281E+09 
0. 2302E+09 
0. 2321E+09 
0. 2339E+09 
0. 2354E+09 
0. 2367E+09 
0 . 2378E+09 
0. 2387E+09 
0 . 2394E+09 
0 . 2399E+09 
0 . 2402E+09 
0. 2404E+09 
0 . 2404E+09 
0. 2402E+09 
0.2 399E+09 
0 . 2395E+09 


0.7513E+08 
0.7966E+08 
0.8440E+08 
0.8936E+08 
0.9453E+08 
0 . 9991E+08 
0 . 1055E+09 
0 . 1113E+09 
0.1173E+09 
0 . 1235E+09 
0 . 1299E+09 
0.1365E+09 
0 . 1432E+09 
0.1501E+09 
O.V^SE+Og 
u .  it>  j  /E+09 
0 . 1709E+09 
0 . 1783E+09 
0.1858E+09 
0 . 1933E+09 
0.2009E+09 
0.2085E+09 
0.2161E+09 
0 . 2236E+09 
0.2311E+09 
0 . 2385E+09 
0.2457E+09 
0.2529E+09 
0.2599E+09 
0.2667E+09 
0.2733E+09 
0.2798E+09 
0.2859E+09 
0.2918E+09 
0.2975E+09 
0.3028E+09 
0. 3079E+09 
0. 3127E+09 
0.3171E+09 
0.3212E+09 
0 . 3250E+09 
0.3285E+09 
0.3317E+09 
0 . 3345E+09 
0.3370E+09 
0.3392E+09 
0.3411E+09 
0 . 3427E+09 
0 . 3440E+09 
0.3450E+09 
0 . 3457E+09 
0 . 3462E+09 
0 . 3464E+09 
0 . 3463E+09 
0.3460E+09 
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0.7050E-02 
0.7100E-02 
0.7150E-02 
0.7200E-02 
0.7250E-02 
0.7300E-02 
0.7350E-02 
0.7400E-02 
0.7450E-02 
0.7500E-02 
0.7550E-02 
0.7600E-02 
0.7650E-02 
0.7700E-02 
0.7750E-02 
0. 7800E-02 
0.7850E--02 
0.7900E-02 
0.7950E-02 
0.8000E-02 
0.8050E-02 
0.8100E-02 
0.8150E-02 
0.8200E-02 
0.8250E-02 
0. 8300E-02 
0.8350E-02 
0.8400E— 02 
0.8450E-02 
0.8500E-02 
0.8550E-02 
0. 8600E-02 
0. 8650E-02 
0.8700E-02 
0.8750E-02 
0.8800E-02 
0. 8850E-02 
0.8900E-02 
0. 8950E-02 
0.9005E-02 
0.9055E-02 
0.9105E-02 
0.9155E-02 
0. 9205E-02 
0 . 9255E-02 
0 .  9305E-02 
0. 9355E-02 
0.9405E-02 
0 . 9455E-02 
0. 9505E-02 
0.9555E-02 
0.9605E-02 
0.9655E-02 
0.9705E-02 
0 . 9755E-02 


0. 3089E+06 
0. 3080E+06 
0. 3069E+06 
0.3057E+06 
0. 3043E+06 
0.3029E+06 
0. 3013E+06 
0. 2996E+06 
0 . 2978E+06 
0 . 2959E+06 
0.2939E+06 
0.2919E+06 
0. 2898E+06 
0.2876E+06 
0 . 2854E+06 
0 . 2831E+06 
0.2808E+06 
0.2785E+06 
0. 2761E+06 
0. 2737E+06 
0.2713E+06 
0.2688E+06 
0.2664E+06 
0.2639E+06 
0.2614E+06 
0 . 259QE+06 
0 . 2563E+06 
0.2527E+06 
0 , 2488E+06 
0.2447E+06 
0. 2404E+06 
0. 2361E+06 
0.2316E+06 
0. 2271E+06 
0.2226E+06 
0 . 2181E+06 
0. 2137E+06 
0.2093E+06 
0. 2049E+06 
0 . 2002E+06 
0 . 1959E+06 
0. 1918E+06 
0 . 1877E+06 
0 . 1837E+06 
0. 1798E+06 
0. 1760E+06 
0 . 1722E+06 
0. 1686E+06 
0 . 1651E+06 
0 . 1617E+06 
0 . 1583E+06 
0. 1S51E+06 
0 . 1519E+06 
0. 1488E+06 
0. 1458E+06 


0.6753E+03 
0.6907E+03 
0.7061E+03 
0.7214E+03 
0.7367E+03 
0.7518E+03 
0.7669E+03 
0.7820E+03 
0.7969E+03 
0.8117E+03 
0.8265E+03 
0.8411E+03 
0.8557E+03 
0.8701E+03 
0.8844E+03 
0.8986E+03 
0.9127E+03 
0.9267E+03 
0.9406E+03 
0.9543E+03 
0.9680E+03 
0.9815E+03 
0.9948E+03 
0. 1008E+04 
0. 1021E+04 
0. 1034E+04 
0. 1047E+04 
0. 1060E+04 
0. 1072E+04 
0. 1085E+04 
0. 1097E+04 
0.1109E+04 
0.1120E+04 
0.1132E+04 
0. 1143E+04 
0. 1154E+04 
0. 1165E+04 
0. 1176E+04 
0, 1186E+04 
0. 1197E+04 
0. 1207E+04 
0. 1217E+04 
0. 1226E+04 
0. 1235E+04 
0. 1245E+04 
0. 1253E+04 
0. 1262E+04 
0. 1271E+04 
0. 1279E+04 
0 . 1287E+04 
0. 1295E+04 
0.1303E+04 
0. 1311E+04 
0. 1318E+04 
0. 1326E+04 


0.8863E+00 
0.9205E+00 
0.9554E+00 
0.9911E+00 
0.1028E+01 
0.1065E+01 
0. 1103E+01 
0.1141E+01 
0. 1181E+01 
0. 1221E+01 
0.1262E+01 
0. 1304E+01 
0. 1346E+01 
0. 1389E+01 
0. 1433E+01 
0. 1478E+01 
0. 1523E+01 
0. 1569E+01 
0. 1616E+03L 
0. 1663E+01 
0. 1711E+01 
0. 1760E+01 
0. 1809E+01 
0. 1859E+01 
0.1910E+01 
0.1962E+01 
0. 2014E+01 
0.2066E+01 
0.2120E+01 
0.2173E+01 
0.2228E+01 
0. 2283E+01 
0.2339E+01 
0. 2395E+01 
0.2452E+01 
0. 2509E+01 
0. 2567E+01 
0.2626E+01 
0.2685E+01 
0. 2751E+01 
0. 2811E+01 
0. 2871E+01 
0. 2932E+01 
0.2994E+01 
0. 3056E+01 
0. 3118E+01 
0 . 3181E+01 
0.3245E+01 
0 . 3308E+01 
0. 3372E+01 
0, 3437E+01 
0 . 3502E+01 
0. 3567E+01 
0. 3633E+01 
0. 3699E+01 


0.3120E+09 
0. 3111E+09 
0.3101E+09 
0. 3089E+09 
0.3075E+09 
0.3060E+09 
0.3045E+09 
0. 3027E+09 
0. 3009E+09 
0. 2990E+09 
0. 2970E+09 
0.2950E+09 
0.2928E+09 
0.2906E+09 
0.2884E+09 
0. 2860E+09 
0.2837E+09 
0.2813E+09 
0 . 2789E+09 
0.2764E+09 
0.2740E+09 
0.2715E+09 
0.2690E+09 
0. 2664E+09 
0. 2639E+09 
0.2614E+09 
0.2586E+09 
0.2551E+09 
0.2511E+09 
0.2470E+09 
0.2426E+09 
0.2382E+09 
0. 2337E+09 
0. 2292E+09 
0. 2246E+09 
0.2201E+09 
0. 2156E+09 
0 . 2111E+09 
0. 2067E+09 
0 . 2019E+09 
0. 1976E+09 
0.1934E+09 
0. 1893E+09 
0. 1853E+09 
0. 1813E+09 
0. 1774E+09 
0. 1737E+09 
0. 1700E+09 
0.1665E+09 
0. 1630E+09 
0. 1596E+09 
0 . 1563E+09 
0. 1531E+09 
0. 1500E+09 
0. 1470E+09 


0. 2389E+09 
0. 2382E+09 
0. 2373E+09 
0. 2364E+09 
0. 2353E+09 
0. 2342E+09 
0. 2330E+09 
0. 2317E+09 
0. 2303E+09 
0. 2288E+09 
0. 2273E+09 
0.2257E+09 
0. 2241E+09 
0. 2224E+09 
0. 2207E+09 
0.2189E+09 
0. 2172E+09 
0.2153E+09 
0. 2135E+09 
0. 2117E+09 
0.2098E+09 
0. 2079E+09 
0. 2060E+09 
0. 2041E+09 
0. 2022E+09 
0. 2003E+09 
0.1982E+09 
0. 1954E+09 
0. 1924E+09 
0. 1892E+09 
0. 1859E+09 
0.1825E+09 
0.1791E+09 
0. 1756E+09 
0. 1722E+09 
0. 1687E+09 
0. 1652E+09 
0. 1618E+09 
0. 1585E+09 
0. 1548E+09 
0. 1515E+09 
0. 1483E+09 
0. 1451E+09 
0. 1421E+09 
0. 1390E+09 
0. 1361E+09 
0. 1332E+Q9 
0. 1304E+09 
0. 1277E+09 
0 . 1250E+09 
0 . 1224E+09 
0. 1199E+09 
0. 1175E+09 
0 . 1151E+09 
0. 1128E+09 


0 . 3455E+09 
0.3448E+09 
0 . 3439E+09 
0. 3429E+09 
Q.3416EX09 
0.3402E+09 
0.3386E+09 
0.3369E+09 
0.3351E+09 
0.3331E+09 
0.3310E+09 
0.3289E+09 
0. 3266E+09 
0.3243E+09 
0. 3219E+09 
0. 3194E+09 
0.3169E+09 
0. 3143E+09 
0. 3117E+09 
0. 3090E+09 
0.3063E+09 
0. 3036E+09 
0. 3009E+09 
0 . 2981E+09 
0.2954E+09 
0 . 2926E+09 
0. 2895E+09 
0.2856E+09 
0.2812E+09 
0.2766E+09 
0.2718E+09 
0.2669E+09 
0.2619E+09 
0.2568E+09 
0. 2518E+09 
0.2467E+09 
0.2417E+09 
0.2367E+09 
0. 2318E+09 
0 . 2265E+09 
0. 2217E+09 
0. 2170E+09 
0. 2124E+09 
0.2073E+09 
0 . 2034E+09 
0.1991E409 
0 . 1949E+09 
0.1908E+09 
0. 1868E+09 
0. 182SE+09 
0 . 1792E+09 
0 . 1755E+09 
0 . 1719E+09 
0.1684E+09 
0 . 1650E+09 
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0.9805E-02  0. 1429E+06  0.1333E+04  0.3766E+01  0.1441E+09  0.1105E+Q9  0.1617E+0Q 

0.9855E-02  0.1401E+06  0.1340E+04  0.3832E+01  0.1412E+09  0.1083E+09  0.1585E+09 

0.9905E-02  0. 1373E+06  0.1347E+04  0.3900E+01  0.1384E+09  0.1Q62E+09  0.1554E+09 

0.9955E-02  0.1346E+06  0.1354E+04  0.3967E+01  0.1357E+09  G.1041E+09  0.1523E+09 

O.iOOQE-Ol  0.1320E+06  0.1360E+04  0.4035E+01  0.1331E+09  0.1021E+09  0.1494E+09 

0.1005E-01  0. 1295E+06  0.1367E+04  0.4103E+01  0.1305E+09  0.1001E+09  0.1465E+09 

0.1010E-01  0. 1270E+06  0.1373E+04  0.4172E+01  0.1280E+09  0.9821E+08  0.1437E+09 

0 . 1015E-01  0.1246E+06  0.1379E+04  0.4240E+01  0.1256E+09  0.9635E+08  0.1409E+09 

0.1020E-01  0. 1222E+06  0.1386E+04  0.4309E+01  0.1232E+09  0.9453E+08  0.1383E+09 

0.1025E-01  0. 1200E+06  0.1392E+04  0.4379E+01  0.1209E+09  0.9277E+08  0.1357E+09 

0.1030E-01  0. 1177E+06  0.1398E+04  0.4449E+01  0.1186E+09  0.9105E+08  0.1332E+09 

0 . 1035E-01  0.1156E+06  0.1403E+04  0.4519E+01  0.1165E+09  0.8937E+08  0.1307E+09 

deltat  t  0.103950E-01  intg  t  0.103948E-01 
PMAXMEAN  Pa  0. 313603E+09  time  at  FMAXMEAN  S6C  0.685009E-02 

FMAXBASE  Pa  0.240397E+09  time  at  PMAXBASE  sec  0. 682509E-02 

PMAXBREECH  Pa  0.346365E+09  time  at  PMAXBREECH  Sec  0.691510E-02 
muzzle  VELOCITY  m/s  0.140782E+04  time  of  muzzle  velocity  sec  0.103927E-01 
total  initial  energy  available  J  =  0.437119E+08 

FOR  PROPELLANT  1  MASSFRACT  BURNT  IS  0.997470E+00 
total  energy  remaining  in  gas  J=  0.295069E+08 
energy  loss  from  projectile  translation  J-  0.971088E+07 
energy  loss  from  projectile  rotation  J=  O.OOOOOOE+OO 
energy  lost  to  gas  and  propellant  motion  J=  0.292176E+07 
energy  lost  to  bore  resistance  J«*  O.OOOOOOE+OO 
energy  lost  to  recoil  J=  0.182576E-11 
energy  loss  from  heat  transfer  J=  0.146181E+07 
energy  lost  to  air  resistance  J*  O.OOOOOOE+OO 
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